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0 Image superposition processing method. 

0 A plurality of image signals are obtained, which represent a radiation image of a single object or radiation 
images of the single object, and which have different frequency characteristics, each of the image signals being 
made up of a series of image signal components. The image signal components of the plurality of the image 
signals, which image signal components represent corresponding picture elements, are weighted with predeter- 
mined weight factors. The weighted image signal components of the plurality of the image signals, which image 
signal components represent corresponding picture elements, are added to one another, and an addition signal 
is thereby obtained. The value of the weight factor with respect to the frequency components, which have a low 
signal-to-noise ratio, is rendered smaller than the value of the weight factor with respect to the frequency 
components, which have a high signal-to-noise ratio, in accordance with the frequency characteristics of each of 
the image signals. 
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BACKGROUND OF THE INVENTION 
Reld of the Invention 

5 This invention relates to a superposition processing method for a radiation image, wherein an addition 
process is carried out on a plurality of Image signals, which represent a radiation image of a single object or 
radiation images of the single object. This invention also relates to an energy subtraction processing 
method, wherein a subtraction process is canied out on a plurality of image signals representing radiation 
images of a single object. 

10 

Description of the Prior Art 

Techniques for reading out a recorded radiation image in order to obtain an image signal, carrying out 
appropriate image processing on the image signal, and then reproducing a visible image by use of the 

75 processed image signal have heretofore been known in various fields. For example, an X-ray image Is 
recorded on an X-ray film having a small gamma value chosen according to the type of image processing 
to be carried out, the X-ray image is read out from the X-ray film and converted into an electric signal (i.e., 
an image signal), and the image signal is processed and then used for reproducing the X-ray image as a 
visible image on a photocopy, or the like. In this manner, a visible image having good image quality with 

20 high contrast, high sharpness, high graininess, or the like, can be reproduced. 

Further, it has been proposed to use stimulable phosphors in radiation image recording and reproduc- 
ing systems. Specifically, a radiation image of an object, such as a human body, Is recorded on a sheet 
provided with a layer of the stimulable phosphor (hereinafter referred to as a stimulable phosphor sheet). 
The stimulable phosphor sheet, on which the radiation image has been stored, is then scanned with 

25 stimulating rays, such as a laser beam, which cause it to emit light in proportion to the amount of energy 
stored thereon during its exposure to the radiation. The light emitted by the stimulable phosphor sheet, 
upon stimulation thereof, is photoelectrically detected and converted into an electric image signal. The 
image signal is then processed and used for the reproduction of the radiation Image of the object as a 
visible image on a recording material. 

30 Techniques for carrying out superposition processing on radiation images have heretofore been 
disclosed In. for example, U.S. Patent No. 4,356.398. In general, radiation Images are used for diagnoses of 
illnesses and for other purposes. When a radiation image is used for such purposes, it is required that even 
small differences in the radiation energy absorption characteristics among structures of an object can be 
detected accurately In the radiation Image. The extent, to which such differences in the radiation energy 

35 absorption characteristics can be detected in a radiation image, is referred to as the contrast detection 
performance or simply as the detection performance. A radiation Image having better detection performance 
has better image quality and can serve as a more effective tool In, particularly, the efficient and accurate 
diagnosis of an illness. Therefore, in order for the Image quality to be improved, it is desirable that the 
detection performance of the radiation image may be enhanced. Practically, the detection performance is 

40 adversely affected by various noises. 

For example, in radiation image recording and reproducing systems using stimulable phosphor sheets, 
it has been found that the noises described below occur during the step for recording a radiation image on 
a stimulable phosphor sheet and reading out the radiation Image therefrom. 
(1) A quantum noise of radiation produced by a radiation source. 

45 (2) A noise due to nonuniformity in how a stimulable phosphor coated on the stimulable phosphor sheet 
is distributed or how stimulable phosphor grains are distributed on the stimulable phosphor sheet. 

(3) A noise of stimulating rays, which cause the stimulable phosphor sheet to emit light in proportion to 
the amount of energy stored thereon during its exposure to radiation. 

(4) A noise of light, which is emitted by the stimulable phosphor sheet, guided and detected. 
50 (5) An electric noise in the system for amplifying and processing an electric signal. 

Superposition processing Is carried out In order to reduce the aforesaid noises markedly so that even 
small differences in the radiation energy absorption characteristics among structures of an object can be 
found accurately in a visible radiation image, which is reproduced finally, i.e. the detection performance of 
the radiation image can be improved markedly. Ordinary techniques and effects of the superposition 
55 processing are as described below. 

A radiation image is stored on each of a plurality of recording media, which have been placed one upon 
another. Thereafter, an image read-out operation Is carried out for each of the recording media. A plurality 
of image signals, which have been obtained from the image read-out operations, are then superposed one 
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upon another. In this manner, various noises described above can be reduced. Specifically, in general, 
noises described in (1) through (5) for the stimulable phosphor sheets exhibit different distributions for 
different radiation images stored on the stimulable phosphor sheets. When the image signals detected from 
the stimulable phosphor sheets are superposed one upon another, the noises can be averaged. Therefore. 

5 the noises become imperceptible in a superposition image, which Is obtained from superposition process- 
ing. Specifically, an image signal having a high signal-to-noise ratio (S/N ratio) is obtained from super- 
position processing. The same effects can be obtained also when radiation images having been recorded 
on sheets of X-ray film are read out. More specifically, most of the noises described in (1) through (5), 
particularly, the noise described in (1), which is one of dominant factors among the noises in a radiation 

70 image, can be approximated by the Poisson statistics. In cases where noises can be approximated by the 
Poisson statistics and two radiation images yield equivalent levels of signals SI and S2 and equivalent 
levels of noises N1 and N2, the level of the signal corresponding to a superposition image, which is 
obtained by carrying out superposition processing on the two radiation images, becomes equal to S1 + 32. 
and the level of noise in the superposition image is represented by Formula (1). 

75 

VNi2 + (1) 

As for the signal-to-noise ratio, which is one of indexes representing the detection performance of a 
radiation image, the signal-to-noise ratios of the two radiation images prior to superposition processing are 
20 represented by the formulas S1/N1 and S2/N2. After superposition processing has been carried out on the 
two radiation images, the signal-to-noise ratio of the resulting superposition image is represented by 
Formula (2). 



(Si + S2)/VNi2 + N22 (2) 

25 

Therefore, as a result of superposition processing, the signal-to-noise ratio can be improved. When 
superposition processing is carried out on image signals representing the two radiation images, the values 
of the image signals may be weighted such that a markedly high signal-to-noise ratio can be obtained. 

By way of example, when superposition processing is to be carried out by using the stimulable 

30 phosphor sheets, two stimulable phosphor sheets have heretofore been housed in a cassette such that they 
overlap one upon the other. Radiation images of an object are then recorded on the two stimulable 
phosphor sheets housed in the cassette. Thereafter, an image read-out operation is carried out on each of 
the two stimulable phosphor sheets, and two image signals are thereby obtained. 

Also, techniques for carrying out subtraction processing on radiation images have heretofore been 

35 known. When subtraction processing is to be carried out, a plurality of (basically, two) radiation images 
recorded under different conditions are photoelectrically read out, and digital image signals which represent 
the radiation images are thereby obtained. The image signal components of the digital image signals, which 
represent corresponding picture elements in the radiation images, are then subtracted from each other, and 
a difference signal is thereby obtained which represents the image of a specific stnjcture or part of the 

40 object represented by the radiation images. With the subtraction processing method, the plurality of digital 
image signals are subtracted from each other in order to obtain a difference signal, and the radiation image 
of a specific structure can be reproduced from the difference signal. 

Basically, subtraction processing is carried out with either the so-called temporal (time difference) 
subtraction processing method or the so-called energy subtraction processing method. In the former 

45 method, in order to extract the image of a specific structure of an object from the image of the entire object, 
the image signal representing a radiation image obtained without injection of contrast media is subtracted 
from the image signal representing a radiation image in which the image of the specific structure of the 
object is enhanced by the injection of contrast media. In the latter method, an object is exposed to several 
kinds of radiation having different energy distributions. Alternatively, the energy distribution of the radiation 

50 carrying image information of an object, is changed after it has been irradiated onto one of a plurality of 
radiation image recording media, after which the radiation impinges upon the second radiation image 
recording medium. In this manner, a plurality of radiation images, in which different images of a specific 
structure of the object are embedded, are obtained. Thereafter, the image signals representing the plurality 
of radiation images are weighted appropriately, when necessary, and subjected to a subtraction process. 

55 and the image of the specific structure of the object is thereby extracted. 

In the aforesaid radiation image recording and reproducing systems utilizing the stimulable phosphor 
sheets, the radiation image stored on the stimulable phosphor sheet is read out directly as an electric 
image signal. Therefore, with such radiation image recording and reproducing systems, the aforesaid 
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subtraction processing can be carried out easily. In cases where energy subtraction processing is to be 
carried out by using the stimulable phosphor sheets, radiation images may be stored on, for example, two 
stimulable phosphor sheets such that the parts of the radiation Images con^esponding to a specific structure 
may be different in the two radiation images. For this purposes, two-shot energy subtraction processing 

6 may be employed wherein the operation for recording a radiation image is carried out twice with two kinds 
of radiation having different energy distributions. Alternatively, one-shot energy subtraction processing may 
be employed wherein, for example, two stimulable phosphor sheets placed one upon the other (they may 
be in contact with each other or spaced away from each other) are simultaneously exposed to radiation, 
which has passed through an object, such that they may be exposed to radiation having different energy 

70 distributions. 

As a method for photoelectrically detecting light emitted by a stimulable phosphor sheet, a method for 
detecting light emitted by two surfaces of a stimulable phosphor sheet has been proposed in, for example, 
U.S. Patent No. 4,346,295. With the proposed method for detecting light emitted by two surfaces of a 
stimulable phosphor sheet, two photoelectric read-out means are located on opposite sides of the 

75 stimulable phosphor sheet. The two surfaces or only one surface of the stimulable phosphor sheet is 
scanned with the stimulating rays, and the light emitted by the two surfaces of the stimulable phosphor 
sheet is photoelectrically detected by the two photoelectric read-out means. With the proposed method for 
detecting light emitted by two surfaces of a stimulable phosphor sheet, a single radiation image is stored on 
the stimulable phosphor sheet, and the light emitted by two surfaces of the stimulable phosphor sheet is 

20 detected on the two sides of the stimulable phosphor sheet. Therefore, the efficiency, with which the light 
emitted by the stimulable phosphor sheet is guided and detected, can be kept high, and a high slgnal-to- 
noise ratio can be obtained. 

With the method for detecting light emitted by two surfaces of a stimulable phosphor sheet, which has 
been proposed in U.S. Patent No. 4.346,295, the stimulable phosphor sheet is placed on a transparent 

25 holder, and two photoelectric read-out means are respectively located above and below the holder. 
Specifically, the light emitted from the front surface of the stimulable phosphor sheet is detected by the 
photoelectric read-out means, which is located above the holder. Also, the light emitted from the back 
surface of the stimulable phosphor sheet is detected by the photoelectric read-out means, which is located 
below the holder. 

30 When the image signals to be subjected to the superposition processing are to be obtained, for 
example, it Is necessary to record radiation images on a plurality of stimulable phosphor sheets superposed 
one upon another. In such cases, the image signal, which is obtained from a stimulable phosphor sheet 
located at the position remote from the radiation source, contains image Information in the low frequency 
band as in the image signal, which is obtained from a stimulable phosphor sheet located at the position 

35 close to the radiation source. However, in the image signal, which is obtained from the stimulable phosphor 
sheet located at the position remote from the radiation source, the frequency dependency in the high 
frequency band is lower than in the image signal, which is obtained from the stimulable phosphor sheet 
located at the position close to the radiation source. In the image signal, which is obtained from the 
stimulable phosphor sheet located at the position remote from the radiation source, as for the high 

40 frequency band, the amount of image information becomes small, and the amount of the noise components 
due to the effects of scattered radiation, or the like, becomes large. Therefore, if the image signal, which is 
obtained from the stimulable phosphor sheet located at the position remote from the radiation source, and 
the image signal, which is obtained from the stimulable phosphor sheet located at the position close to the 
radiation source, are weighted in the same manner and added to each other, the image quality can be kept 

45 good in the low frequency band in the addition signal obtained from the addition process, but the noise 
components will be emphasized and adversely affect the image quality in the high frequency band. Such 
adverse effects upon the image quality will also occur with the image signal, which Is obtained from the 
front surface of the stimulable phosphor sheet by the method for detecting light emitted by two surfaces of 
the stimulable phosphor sheet, and the Image signal, which is obtained from the back surface of the 

50 stimulable phosphor sheet. Further, as for the image signals to be subjected to the energy subtraction 
processing, the proportion of the noise components varies for different frequency bands of the Image signal. 
Therefore, when the subtraction process is carried out on the image signals, it often occurs that the amount 
of the noise components in the difference signal becomes large, depending on the weight factor employed 
for each of the image signals. 

55 Recently, it is desired to increase the speed, with which a radiation image is read out. Therefore, in 
Japanese Unexamined Patent Publication Nos. 60(1985>-117212 and 62{1987)-90615, the applicant pro- 
posed radiation image read-out apparatuses, which are capable of quickly reading out a radiation image 
from a stimulable phosphor sheet. 
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Also, in order for the speed, with which a radiation image is read out, to be increased, it has been 
proposed to increase the response speed of a stimulable phosphor with respect to stimulating rays, I.e. to 
increase the speed of light emission response of the stimulable phosphor with respect to the in'adiation of 
the stimulating rays, with a method wherein, for example, cerium is added to the stimulable phosphor. 

5 When the response speed of the stimulable phosphor sheet with respect to the stimulating rays Is kept 
higher, it becomes necessary to employ stimulating rays having a higher output power for the scanning of 
the stimulable phosphor sheet. Therefore, the stimulable phosphor sheet is quickly scanned with a laser 
beam, which serves as the stimulating rays and which has a high output power of at least 50mW. 

However, in cases where the aforesaid operation for quickly reading out a radiation image is carried out, 

70 the stimulating rays are moved very quickly on the stimulable phosphor sheet. Therefore, with certain kinds 
of stimulable phosphors constituting the stimulable phosphor sheets, the problems often occur in that the 
stimulable phosphor sheet cannot emit light immediately after being exposed to the stimulating rays, and a 
time lag occurs between when the stimulable phosphor sheet is exposed to the stimulating rays and when 
the stimulable phosphor sheet emits light. At an image contour portion, or the like, the amount of the light 

15 emitted by the stimulable phosphor sheet changes sharply. Therefore, the values of the image signal, which 
correspond to an image contour portion, or the like, should change sharply in the main scanning direction. 
However, if the time lag occurs between when the stimulable phosphor sheet is exposed to the stimulating 
rays and when the stimulable phosphor sheet emits light, the values of the obtained Image signal, which 
correspond to an image contour portion, or the like, will not change sharply in the main scanning direction, 

20 As a result, the visible image reproduced from the image signal becomes unsharp in the main scanning 
direction, and the sharpness of the reproduced image cannot be kept high. 

As for an image signal actually obtained from a stimulable phosphor sheet, or the like, the response 
characteristics will vary for different frequency bands. Therefore, such that a well-balanced visible image 
can be reproduced, it is desired to enhance the signal-to-noise ratio by reducing the noise components in 

25 the obtained addition signal or the obtained subtraction signal, and to emphasize or reduce the components 
of a desired frequency band in the obtained addition signal or the obtained subtraction signal. In such 
cases, it is necessary to carry out filtering processing on the obtained addition signal or the obtained 
subtraction signal by using, for example, a mask filter shown in Rgure 35. However, in such cases, the 
problems occur in that a long calculation time is required to carry out the filtering processing, and in that a 

30 complicated apparatus must be used. 

Noise contained in an image signal Is also affected by the dose of radiation delivered to the stimulable 
phosphor sheet. Specifically, when the dose of radiation delivered to the stimulable phosphor sheet is large, 
the proportion of the fixed noise due to the structure of the stimulable phosphor sheet, such as the state in 
which the stimulable phosphor is applied to the stimulable phosphor sheet, becomes higher than the 

35 proportion of the quantum noise of the radiation. Therefore, the ratio, in which a frequency band of the 
image signal is to be weighted such that the image quality of the reproduced image obtained from the 
addition signal or the subtraction signal may be kept best, varies for different doses of radiation delivered to 
the object. 

40 SUMMARY OF THE INVENTION 

The primary object of the present invention is to provide a superposition processing method for a 
radiation image, wherein a superposition image having good image quality and containing little noise 
component is obtained. 

45 Another object of the present invention is to provide an energy subtraction processing method, wherein 
an energy subtraction image having good image quality and containing little noise component is obtained. 

A further object of the present invention is to provide a superposition processing method for a radiation 
image, wherein the addition process on image signals is canried out simply, at a low cost, quickly, and for 
each of different frequencies. 
50 A still further object of the present invention is to provide a superposition processing method for a 
radiation image, wherein a reproduced image having a high sharpness is obtained from image signals 
detected with operations for quickly reading out a radiation image. 

Another object of the present invention is to provide a superposition processing method for a radiation 
image, wherein the addition process on image signals is canried out simply, at a low cost, quickly, and for 
55 each of different frequencies, and wherein the components of a desired frequency band in an image signal 
are capable of being altered. 

A further object of the present invention is to provide an energy subtraction processing method, wherein 
the subtraction process on image signals is carried out simply, at a low cost, quickly, and for each of 
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different frequencies, and wtierein the components of a desired frequency band in an image signal are 
capable of being altered. 

A still further object of the present invention is to provide a superposition processing method for a 
radiation image, wherein a reproduced image having good image quality is obtained regardless of the dose 
s of radiation delivered to an object, and wherein the addition process on Image signals is carried out simply, 
at a low cost, quickly, and for each of different frequencies. 

Another object of the present invention is to provide an energy subtraction processing method, wherein 
a reproduced image having good image quality is obtained regardless of the dose of radiation delivered to 
an object, and wherein the subtraction process on image signals is carried out simply, at a low cost, 
70 quickly, and for each of different frequencies. 

The present invention provides a first superposition processing method for a radiation image, compris- 
ing the steps of: 

i) obtaining a plurality of image signals, which represent a radiation image of a single object or radiation 
images of the single object, and which have different frequency characteristics, each of the image 

75 signals being made up of a series of image signal components, 

ii) weighting the image signal components of the plurality of the image signals, which image signal 
components represent corresponding picture elements, with predetermined weight factors, and 

iii) adding the weighted image signal components of the plurality of the image signals to one another, 
which image signal components represent corresponding picture elements, an addition signal being 

20 thereby obtained, 

wherein the value of the weight factor with respect to the frequency components, which have a low 
signal-to-noise ratio, is rendered smaller than the value of the weight factor with respect to the frequency 
components, which have a high signal-to-noise ratio, in accordance with the frequency characteristics of 
each of the image signals. 

25 The present invention also provides a second superposition processing method for a radiation image, 
wherein the first superposition processing method for a radiation image in accordance with the present 
invention is modified such that the image superposition processing may be carried out by: 

a) carrying out a Fourier transform on each of the image signals, each of the image signals being 
thereby decomposed into a plurality of Fourier transform factor signals, each being of one of different 

30 frequency bands, 

b) weighting the Fourier transform factor signals, which are of a single frequency band, with the weight 
factors, such that the values of the weight factors may be varied for the Fourier transform factor signals 
of the different frequency bands, 

c) adding the weighted Fourier transform factor signals, which are of a single frequency band, to one 
35 another, an addition Fourier transfomn factor signal being thereby obtained for each of the different 

frequency bands, and 

d) carrying out an inverse Fourier transform on the addition Fourier transform factor signal. 

The present invention further provides a third superposition processing method for a radiation image, 
wherein the first superposition processing method for a radiation image in accordance with the present 
40 invention is modified such that the image superposition processing may be carried out by: 

a) subjecting each of the image signals to a transform to multi-resolution space, each of the image 
signals being thereby decomposed into a plurality of transform factor signals, each being of one of 
different frequency bands. 

b) weighting the transform factor signals, which are of a single frequency band, with the weight factors, 
45 such that the values of the weight factors may be varied for the transform factor signals of the different 

frequency bands, 

c) adding the weighted transfoon factor signals, which are of a single frequency band, to one another, an 
addition transform factor signal being thereby obtained for each of the different frequency bands, and 

d) carrying out an inverse transform on the addition transform factor signal. 

50 The term "transfonm to multi-resolution space" as used herein means a transform, such as a wavelet 
transform or a sub-band transform, wherein an image signal is decomposed into a plurality of signals of 
different frequency bands by using a filter, which is shorter than a filter employed in the Fourier transform. 
How the wavelet transform is canried out will be described hereinbelow. 

The wavelet transform has recently been developed as a frequency analysis method and has heretofore 
55 been applied to stereo pattern matching, signal compression, and the like. The wavelet transform is 
described in, for example, "Wavelets and Signal Processing," by Olivier Rioul and Martin Vetterii, IEEE SP 
Magazine, pp. 14-38, October 1991; and "Zero-Crossings of a Wavelet Transform," by Stephane Mallat, 
IEEE Transactions on Information Theory, Vol. 37, No. 4, pp. 1019-1033. July 1991. 
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With the wavelet transform, a signal is transformed into frequency signals, each being of one of a 
plurality of different frequency bands, by utilizing a function h, which is shown in Rgure 22, as a basic 
function and in accordance with the formula 



W (a. b) = 



10 



OO 



f ( t ) h ( a . b ) d t 



(3) 



75 



wherein 

f(t): the signal having an arbitrary wave form, 

W(a,b): the wavelet transform of f(t), 

h(a. b) = vJh(at-b) 



a: the degree of contraction of the function, 

b: the amount of movement in the horizontal axis direction. 

20 Therefore, the problems with regard to a false oscillation, which occurs with the Fourier transfomn, do not 
occur. Specifically, when filtering processing is canied out by changing the period and the degree of 
contraction of the function h and moving the function h on an original signal, frequency signals, each of 
which is adapted to one of desired frequencies ranging from a fine frequency to a coarse frequency. By 
way of example, Figure 23 shows signals, which are obtained by carrying out the wavelet transform on an 

25 original signal Sorg and then can7ing out the inverse wavelet transform for each of frequency bands. Figure 
24 shows signals, which are obtained by carrying out the Fourier transform on the original signal Sorg and 
then carrying out the inverse Fourier transform for each of the frequency bands. As will be understood from 
Figures 23 and 24, the wavelet transform has the advantage over the Fourier transform in that a frequency 
signal of a frequency band corresponding to the oscillation of the original signal Sorg can be obtained. 

30 Specifically, with the Fourier transform, an oscillation occurs in a part B' of a frequency band 7, which 
corresponds to a part B of the original signal Sorg. However, with the wavelet transform, as in the original 
signal Sorg, no oscillation occurs in a part A* of a frequency band W7, which corresponds to a part A of the 
original signal Sorg, 

The present invention still further provides a fourth superposition processing method for a radiation 
35 image, wherein the third superposition processing method for a radiation image in accordance with the 
present invention is modified such that the transform to multi-resolution space may be carried out with a 
wavelet transform. 

The present invention also provides a fifth superposition processing method for a radiation image, 
wherein the first, third, or fourth superposition processing method for a radiation image in accordance with 

40 the present invention is modified such that the value of the weight factor with respect to a portion of the 
radiation image, to which a large dose of radiation reached during an operation for recording the radiation 
image, may be rendered larger than the value of the weight factor with respect to a portion of the radiation 
image, to which a small dose of radiation reached during the operation for recording the radiation image, in 
accordance with the portions of the object, the patterns of which are embedded in the radiation image. 

45 The present invention further provides a sixth superposition processing method for a radiation image, 
wherein the first, second, third, fourth, or fifth superposition processing method for a radiation image in 
accordance with the present invention is modified such that the plurality of the image signals are obtained 
by: 

exposing the two surfaces or either one of the two surfaces of a single stimulable phosphor sheet, on 
50 which the radiation image has been stored, to stimulating rays, which cause the stimulable phosphor sheet 
to emit light in proportion to the amount of energy stored thereon during its exposure to radiation, and 

photoelectrically detecting the emitted tight independently on the opposite surface sides of the 
stimulable phosphor sheet. 

The present invention still further provides a seventh superposition processing method for a radiation 
55 image, wherein the first superposition processing method for a radiation image in accordance with the 
present invention is modified such that the plurality of the image signals may be analog image signals, and 
the image superposition processing may be carried out by: 
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a) filtering alt of the analog inDage signals with filters, which have weights changing the frequency 
characteristics of the analog innage signals, and 

b) adding the analog image signals, which have been otained from the filtering process, to one another. 
The present invention also provides an eighth superposition processing method for a radiation image, 

5 wherein the first superposition processing method for a radiation image in accordance with the present 
invention is modified such that the plurality of the image signals may be analog image signals, and the 
image superposition processing may be carried out by: 

a) filtering at least a single desired analog image signal, which is among the analog image signals, with a 
filter, which has a weight changing the frequency characteristics of the desired analog image signal, and 
10 b) adding the analog image signal, which has been obtained from the filtering process, and the other 
analog image signals to one another. 

The present invention further provides a ninth superposition processing method for a radiation image, 
wherein the first, second, third, fourth, fifth, seventh, or eighth superposition processing method for a 
radiation image in accordance with the present invention is modified such that the plurality of the image 
75 signals are obtained by: 

exposing each of at least two stimulable phosphor sheets, on each of which a radiation image has been 
stored, to stimulating rays, which cause the stimulable phosphor sheet to emit light in proportion to the 
amount of energy stored thereon during its exposure to radiation, and 
photoelectrically detecting the emitted light. 
20 The present invention still further provides a first energy subtraction processing method comprising the 
steps of: 

i) forming a plurality of radiation images of a single object respectively with a plurality of kinds of 
radiation having different energy distributions, different images of at least part of the object being 
embedded in the plurality of the radiation images, 
25 ii) obtaining a plurality of image signals, each of which is made up of a series of image signal 
components, from the plurality of the radiation images, 

iii) weighting the image signal components of the plurality of the image signals, which image signal 
components represent corresponding picture elements in the plurality of the radiation images, with 
predetermined weight factors, and 
30 iv) subtracting the weighted image signal components of the plurality of the image signals from one 
another, which image signal components represent corresponding picture elements in the plurality of the 
radiation images, a difference signal representing an image of a specific structure of the object being 
thereby obtained, 

wherein the value of the weight factor with respect to the frequency components, which have a low 
35 signal-to-noise ratio, is rendered smaller than the value of the weight factor with respect to the frequency 
components, which have a high signal-to-noise ratio, in accordance with the frequency characteristics of 
each of the image signals. 
In the above-described first energy subtraction processing method and also in the below-described 
various other methods according to the present invenfion, including those defined in claims, the expression 
40 of "with a plurality of kinds of radiafion having different energy distributions" does not necessarily means a 
plurality of separate radiations but includes a plurality of kinds of radiations originated from a single 
radiafion wherein, for example, one is a direct radiation fi'om a radiation source, and the other is a radiation 
from the same radiation source, which has passed through a recording medium (e.g., a stimulable phosphor 
sheet) and/or a filter, or the like, and the low energy components of which have been filtered out. Therefore, 
45 the plurality of the radiation images can be formed one after another by using different radiafions having 
different energy distribufions. Alternatively, the plurality of the radiation images can be formed simulta- 
neously by using a single radiation and placing a plurality of recording media (e.g., stimulable phosphor 
sheets) one upon anotiier with or without a filter interposed therebetween. When the fitter is not used, the 
stimulable phosphor sheet located closer to the radiation source serves as a filter for filtering out the tow 
50 energy components of the radiation. 

The present invention also provides a second energy subtraction processing method, wherein the first 
energy subtraction processing method in accordance with the present invention is modified such that the 
difference signal may be obtained by: 

a) carrying out a Fourier transform on each of the image signals, each of the image signals being 
55 thereby decomposed into a plurality of Fourier transform factor signals, each being of one of different 

frequency bands, 

b) weighting the Fourier transform factor signals, which are of a single frequency band, with the weight 
factors, such that the values of the weight factors may be varied for the Fourier transform factor signals 
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of the different frequency bands, 

c) subtracting the weighted Fourier transform factor signals, which are of a single frequency band, from 
one another, a subtraction Fourier transform factor signal being thereby obtained for each of the different 
frequency bands, and 

6 d) carrying out an inverse Fourier transform on the subtraction Fourier transform factor signal. 

The present invention further provides a third energy subtraction processing method, wherein the first 
energy subtraction processing method in accordance with the present invention is modified such that the 
subtraction signal may be obtained by: 

a) subjecting each of the image signals to a transform to multi-resolution space, each of the image 
10 signals being thereby decomposed into a plurality of transform factor signals, each being of one of 

different frequency bands, 

b) weighting the transform factor signals, which are of a single frequency band, with the weight factors, 
such that the values of the weight factors may be varied for the transform factor signals of the different 
frequency bands, 

75 c) subtracting the weighted transform factor signals, which are of a single frequency band, from one 
another, a subtraction transform factor signal being thereby obtained for each of the different frequency 
bands, and 

d) carrying out an inverse transform on the subtraction transform factor signal. 

The present invention still further provides a fourth energy subtraction processing method, wherein the 
20 third energy subtraction processing method in accordance with the present invention is modified such that 
the transform to multi-resolution space may be carried out with a wavelet transform. 

The present invention also provides a fifth energy subtraction processing method, wherein the first, 
third, or fourth energy subtraction processing method in accordance with the present invention is modified 
such that the value of the weight factor with respect to a portion of each radiation image, to which a large 
25 dose of radiation reached during an operation for recording the radiation image, may be rendered larger 
than the value of the weight factor with respect to a portion of the radiation image, to which a small dose of 
radiation reached during the operation for recording the radiation image, in accordance with the portions of 
the object, the patterns of which are embedded in the radiation image. 

The present invention further provides a sixth energy subtraction processing method, wherein the first 
30 energy subtraction processing method in accordance with the present invention is modified such that the 
plurality of the image signals may be analog image signals, and the difference signal may be obtained by; 

a) filtering all of the analog image signals with filters, which have weights changing the frequency 
characteristics of the analog image signals, and 

b) subtracting the analog image signals, which have been obtained from the filtering process, from one 
35 another. 

The present invention still further provides a seventh energy subtraction processing method, wherein 
the first energy subtraction processing method in accordance with the present invention is modified such 
that the plurality of the image signals may be analog image signals, and the difference signal may be 
obtained by: 

40 a) filtering at least a single desired analog image signal, which is among the analog image signals, with a 
filter, which has a weight changing the frequency characteristics of the desired analog image signal, and 
b) subtracting the analog image signal, which has been obtained from the filtering process, and the other 
analog image signals from one another. 

With the first to ninth superposition processing methods for a radiation Image in accordance with the 
45 present invention, the value of the weight factor with respect to the frequency components, which have a 
low signal-to-noise ratio, is rendered smaller than the value of the weight factor with respect to the 
frequency components, which have a high signal-to-noise ratio. In accordance with the frequency char- 
acteristics of each of the image signals. Therefore, the addition signal, which is obtained from the addition 
process, has a high signal-to-noise ratio over the entire frequency bands. Accordingly, a superposition 
50 image having good image quality can be reproduced from the addition signal. 

Specifically, the Fourier transform may be carried out on each of the image signals, and each of the 
image signals may thereby be decomposed into a plurality of Fourier transform factor signals, each being of 
one of different frequency bands. The weighted addition process may then be carried out on the Fourier 
transform factor signals, which are of a single frequency band, such that the values of the weight factors 
55 may be varied for the Fourier transform factor signals of the different frequency bands. In this manner, the 
value of the weight factor with respect to the frequency components, which have a low signal-to- noise ratio, 
can be rendered smaller than the value of the weight factor with respect to the frequency components, 
which have a high signal-to-noise ratio. Therefore, a superposition image having good image quality can be 
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obtained. 

Also, each of the image signals may be subjected to the transform to multi-resolution space, such as 
the wavelet transform or the sub-band transform, and may thereby be decomposed into a plurality of the 
transform factor signals, each being of one of different frequency bands. An addition signal may then be 
5 obtained from the transform factor signals. In this manner, the image signal can be decomposed into the 
plurality of the frequency bands by using a short filter. Therefore, the apparatus for carrying out the 
superposition processing method for a radiation image in accordance wrth the present invention can be kept 
simple. 

Further, the transform factor signals, each of which is of one of different frequency bands, are obtained 

70 by transforming each of the image signals to multi-resolution space. Each of the transform factor signals 
thus obtained is constituted of an image signal contracted from the image signal before being transformed. 
Therefore, the dose of radiation, which reached to a portion of the radiation image during the operation for 
recording the radiation image, can be detected by carrying out. for example, an analysis of the probability 
density function of each of the transform factor signals. Thereafter, the value of the weight factor with 

/5 respect to a portion of the radiation image, to which a large dose of radiation reached during the operation 
for recording the radiation image, can be rendered larger than the value of the weight factor with respect to 
a portion of the radiation image, to which a small dose of radiation reached during the operation for 
recording the radiation image. In this manner, an image having good image quality can be obtained. 

Furthermore, the plurality of the image signals may be taken as analog Image signals. Filtering may be 

20 carried out on all of the analog image signals by using the filters, which have weights changing the 
frequency characteristics of the analog image signals. In this manner, the weighting process is effected on 
desired frequency bands of the analog image signals. Alternatively, filtering may be carried out on at least a 
single desired analog image signal, which is among the analog image signals, by using a filter, which has a 
weight changing the frequency characteristics of the desired analog image signal. In this manner, the 

25 weighting process is effected on a desired frequency band of the analog image signal. Thereafter, the 
addition signal may be obtained by carrying out the addition process on the image signals. In this manner, 
as In the aforesaid method utilizing the Fourier transform, the wavelet transform, or the sub-band transform, 
a superposition image having good image quality can be obtained from the addition signal. 

In cases where the image signals are obtained with the aforesaid method for detecting light emitted by 

30 two surfaces of a stimulable phosphor sheet, the high frequency components of the image signal, which has 
been obtained from the back surface side of the stimulable phosphor sheet (i.e. the side of the stimulable 
phosphor sheet, which was remote from the radiation source during the operation for recording the radiation 
image), contain a large amount of noise components due to scattered radiation, or the like. Also, in cases 
where the image signals are obtained by recording the radiation images respectively on a plurality of 

35 stimulable phosphor sheets placed one upon another, the high frequency components of the image signal 
having been obtained from a stimulable phosphor sheet, which was remote from the radiation source during 
the operation for recording the radiation images, contain a large amount of noise components due to 
scattered radiation, or the like. Therefore, the value of the weight factor with respect to the high frequency 
components of the image signal, which has been obtained from the side of the stimulable phosphor sheet 

40 remote from the radiation source or which has been obtained from the stimulable phosphor sheet remote 
from the radiation source, is rendered smaller than the value of the weight factor with respect to the high 
frequency components of the image signal, which has been obtained from the side of the stimulable 
phosphor sheet close to the radiation source or which has been obtained from the stimulable phosphor 
sheet close to the radiation source. In this manner, a superposition image containing little noise component 

45 can be obtained. 

The aforesaid process for varying the value of the weight factor for different frequency bands of the 
Image signal can be applied to the energy subtraction processing method. Specifically, the value of the 
weight factor with respect to the frequency components, which have a low signal-to-noise ratio, is rendered 
smaller than the value of the weight factor with respect to the frequency components, which have a high 
50 signal-to-noise ratio, in accordance with the frequency characteristics of each of the image signals, which 
are to be subjected to the subtraction process. In such cases, a subtraction signal, which has been obtained 
from the subtraction process, contains little noise component, and a subtraction image having good image 
quality can be reproduced from the subtraction signal. 

The present invention also provides, a tenth superposition processing method for a radiation image, 
55 comprising the steps of: 

i) obtaining a plurality of image signals, which represent a radiation image of a single object or radiation 
images of the single object, and which have different frequency characteristics, each of the image 
signals being made up of a series of image signal components, and 
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ii) adding the image signal components of the plurality of the image signals to one another, which image 
signal components represent corresponding picture elements, an addition signal being thereby obtained. 

wherein, such that the signal-to-noise ratio of the addition signal may be kept high, the addition 
signal is obtained by: 

5 a) subjecting at least a single desired image signal, which is among the plurality of the Image signals, to 
image processing, which changes the frequency characteristics of the desired image signal, and 

b) adding the image signal, which has been obtained from the image processing, and the other image 
signals to one another. 

The image processing may be carried out on all of the plurality of the image signals. 
70 The image processing may be the processing for carrying out the convolution of the desired image 
signal with a mask filter, which has predetermined frequency characteristics. The processing with the 
convolution may be carried out on all of the plurcility of the image signals. 

In cases where the processing with the convolution is carried out on all of the plurality of the image 
signals, the sum of the frequency characteristics of the mask filters employed for the plurality of the image 
15 signals may be equal to 1 at an arbitrary frequency. 

The tenth superposition processing method for a radiation image in accordance with the present 
invention should preferably be modified such that each of the plurality of the image signals may be 
obtained by scanning a stimulable phosphor sheet, on which the radiation image has been stored, with 
stimulating rays, which cause the stimulable phosphor sheet to emit light in proportion to the amount of 
20 energy stored thereon during its exposure to radiation, in the main scanning direction and the sub-scanning 
direction, the emitted light being detected photoelectrical ly, and 

the image processing is the processing for emphasizing the high frequency components of the desired 
image signal with respect to the main scanning direction. 

The processing should more preferably be the processing for carrying out the convolution of the 
25 desired image signal with respect to the main scanning direction with a mask filter, which has frequency 
characteristics capable of emphasizing the high frequency components of the desired image signal with 
respect to the main scanning direction. 

The image processing with the convolution may be carried out on each of the image signals by using a 
single mask filter. 

30 The present invention further provides an eleventh superposition processing method for a radiation 
image, comprising the steps of: 

i) obtaining two image signals, which represent a radiation image of a single object or radiation images of 
the single object, and which have different frequency characteristics, each of the image signals being 
made up of a series of image signal components, and 
35 ii) adding the image signal components of the two image signals to one another, which image signal 
components represent corresponding picture elements, an addition signal being thereby obtained, 

wherein, such that the signal-to-noise ratio of the addition signal may be kept high, the addition 
signal is obtained by: 

a) obtaining a difference signal between the two image signals, 
40 b) carrying out the convolution of the difference signal with either one of two mask filters, which have 
predetermined frequency characteristics, the sum of the frequency characteristics of the two mask filters 
being equal to 1 at an arbitrary frequency, and 

c) adding the difference signal, which has been obtained from the convolution, and either one of the two 
image signals to each other. 

45 The present invention still further provides an eighth energy subtraction processing method comprising 
the steps of: 

1) forming a plurality of radiation images of a single object respectively with a plurality of kinds of 
radiation having different energy distributions, different images of at least part of the object being 
embedded in the plurality of the radiation images, 
50 ii) obtaining a plurality of image signals, each of which is made up of a series of image signal 
components, from the plurality of the radiation images. 

iii) weighting the image signal components of the plurality of the image signals, which Image signal 
components represent corresponding picture elements in the plurality of the radiation images, with 
predetermined weight factors, and 

55 iv) subtracting the weighted image signal components of the plurality of the image signals from one 
another, which image signal components represent corresponding picture elements in the plurality of the 
radiation images, a difference signal representing an image of a specific structure of the object being 
thereby obtained. 
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wherein, such that the signal-to-noise ratio of the difference signal may be kept high, the difference 
signal is obtained by: 

a) subjecting at least a single desired image signal, which is among the plurality of the image signals, to 
image processing, which changes the frequency characteristics of the desired image signal, and 
5 b) subtracting the image signal, which has been obtained from the image processing, and the other 
image signals from one another. 

In the eighth energy subtraction processing method in accordance with the present invention, the image 
processing may be carried out on all of the plurality of the image signals. 

Also, the image processing may be the processing for carrying out the convolution of the desired image 

70 signal with a mask filter, which has predetermined frequency characteristics. The processing with the 
convolution may be carried out on all of the plurality of the image signals. 

With the tenth superposition processing method for a radiation image in accordance with the present 
invention, at least a single desired image signal, which is among the plurality of the image signals, is 
subjected to the image processing, which changes the frequency characteristics of the desired image 

15 signal. The frequency characteristics of the desired image signal can be changed such that noise may be 
reduced. Therefore, the radiation image represented by the addition signal can have good image quality 
and little noise. In addition, with the tenth superposition processing method for a radiation image in 
accordance with the present invention, the image processing, which changes the frequency characteristics 
of the image signal, is carried out on the entire image signal. Therefore, it is not necessary to carry out a 

20 frequency transform, such as the wavelet transform or the Fourier transform. Accordingly, the amount of 
calculation can be kept small, and the apparatus for carrying out the superposition processing method for a 
radiation image in accordance with the present invention can be kept simple. As a result, a superposition 
image having good image quality can be obtained quickly and at a low cost. 

Also, in cases where the image processing, which changes the frequency characteristics, is carried out 

25 on alt of the plurality of the image signals, an addition signal having better image quality can be obtained. 

Further, in cases where the image processing is carried out by the convolution of the image signal with 
a mask filter, addition signals, in which the responses with respect to various frequency bands have been 
emphasized, can be obtained by changing the frequency characteristics of the mask filter. 

Furthermore, in cases where the processing with the convolution is earned out on all of the plurality of 

30 the image signals, an addition signal having better image quality can be obtained. 

Moreover, in cases where the sum of the frequency characteristics of the mask filters, which are 
employed for the plurality of the image signals, is equal to 1 at an arbitrary frequency, when the image 
signals obtained from the image processing are added to one another, it becomes unnecessary for the 
weighting process to be canned out such that the addition ratio of the image signals may be equal to 1. 

35 Therefore, the operation time can be kept short, and the addition process can be carried out quickly. 

Also, in cases where the image processing with the convolution is carried out on each of the image 
signals by using a single mask filter, the number of the mask filters, which are to be stored in an apparatus 
for carrying out the superposition processing method for a radiation image in accordance with the present 
invention, can be kept small. Therefore, the apparatus can be kept simple. 

40 With the eleventh superposition processing method for a radiation image in accordance with the present 
invention, wherein the superposition processing is carried out on two Image signals, the difference signal 
between the two image signals is obtained. The convolution of the difference signal is carried out with either 
one of two mask filters, the sum of the frequency characteristics of the two mask filters being equal to 1 at 
an arbitrary frequency. The difference signal, which has been obtained from the convolution, and either one 

45 of the two image signals are then added to each other. Specifically, the process is carried out which is 
represented by one of the formulas 

Sadd = Si + F2 "(Sa - Si) 
Sadd = S2 + Fi - (Si -S2) 

50 

where S1 and S2 represent the image signals, Sadd represents the addition signal, and F1 and F2 
represent the mask filters satisfying the condition of F1 +F2 = 1. Therefore, only a single convolution may 
be canied out on the values of the difference between the image signals. Also, the number of the mask 
filters, which are to be stored in the apparatus for carrying out the superposition processing method for a 
55 radiation image in accordance with the present invention, can be kept small. Therefore, the operation time 
can be kept short, and the apparatus can be kept simple. 

In the tenth superposition processing method for a radiation image in accordance with the present 
invention, each of the plurality of the Image signals may be obtained by scanning a stimulable phosphor 
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sheet, on which the radiation innage has been stored, with stimulating rays, which cause the stimulable 
phosphor sheet to emit light in proportion to the amount of energy stored thereon during its exposure to 
radiation, in the main scanning direction and the sub-scanning direction, the emitted light being detected 
photoelectrically. In such cases, the image processing may be the processing for carrying out the 

5 convolution of the desired image signal with respect to the main scanning direction with a mask filter, which 
has the frequency characteristics capable of emphasizing the high frequency components of the desired 
image signal with respect to the main scanning direction. Also, as the image processing, the convolution of 
the desired image signal may be carried out with respect to the sut)-scanning direction with a mask filter, 
which has predetermined frequency characteristics. In this manner, the high frequency components of the 

70 image signal with respect to the main scanning direction, i.e. an image part, such as an contour part of the 
image, at which the values of the image signal change sharply, can be emphasized. As a result, the part, at 
which the value of the image signal should change sharply, but which has been rendered unsharp due to 
the quick image read-out operation, can be emphasized. An image signal, in which the problems with 
regard to the unsharp change have been eliminated, can thus be obtained. Therefore, in cases where the 

75 image read-out operation is carried out quickly, the image signal can be prevented from becoming unsharp 
with respect to the main scanning direction, and a reproduced image having a high sharpness can be 
obtained. 

The processing carried out in the tenth or eleventh superposition processing method for a radiation 
image in accordance with the present invention can also be applied to the energy subtraction processing in 
20 the eighth energy subtraction processing method in accordance with the present invention. Therefore, with 
the eighth energy subtraction processing method in accordance with the present invention, a difference 
signal representing an image having good image quality and containing little noise can be obtained from the 
subtraction processing. Also, the image processing, which changes the frequency characteristics of the 
image signal, is carried out on the entire image signal. Therefore, it is not necessary to carry out a 
25 frequency transform, such as the wavelet transform or the Fourier transform. Accordingly, the amount of 
calculation can be kept small, and the apparatus for candying out tiie eighth energy subtraction processing 
method in accordance with ttie present invention can be kept simple. As a result, a subtraction image 
having good image quality can be obtained quickly and at a low cost. 

The present invention still further provides a twelfth superposition processing method for a radiation 
30 image, comprising the steps of: 

i) obtaining a plurality of image signals, which represent a radiation image of a single object or radiation 
images of the single object, and which have different frequency characteristics, each of the image 
signals being made up of a series of image signal components, and 

11) adding the image signal components of the plurality of the image signals to one anotiier, which image 
35 signal components represent corresponding picture elements, an addition signal being thereby obtained, 
wherein the addition signal is obtained by: 

a) carrying out the convolution of at least a single desired image signal, which is among the plurality of 
the image signals, with a mask filter, which has frequency characteristics capable of keeping the signal- 
to-noise ratio of the addition signal high and altering tiie response characteristics of the addition signal 

40 with respect to a desired frequency band when the addition signal is obtained by adding an image 
signal, that Is obtained from the convolution of the desired image signal, and the other image signals to 
one another, a processed image signal being thereby obtained from the convolution, and 

b) adding the processed image signal and the other image signals to one another. 

The twelfth superposition processing method for a radiation image in accordance with the present 
45 invention should preferably be modified such that each of the plurality of the image signals may be 
obtained by scanning a sheet-like recording medium, on which the radiation image has been recorded, with 
a light beam in the main scanning direction and the sub-scanning direction, the radiation image being 
thereby read out photoelectrically, and 

the processed image signal may be obtained by carrying out the convolution of the desired image 
50 signal with the mask filter, which has different frequency characteristics with respect to the main scanning 
direction and the sub-scanning direction, such tiiat the frequency response characteristics of the addition 
signal with respect to the main scanning direction and the frequency response characteristics of the 
addition signal with respect to the sub-scanning direction may become approximately identical with each 
other. 

55 Also, in the twelfth superposition processing method for a radiation image in accordance with the 
present invention, the processing with the convolution may be carried out on all of the plurality of the image 
signals. 
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In cases where the processing with the convolution is carried out on all of the plurality of the image 
signals, the term "other image signals" as used herein also means the image signals obtained from the 
convolution. 

The present invention also provides a ninth energy subtraction processing method comprising the steps 

5 of: 

i) forming a plurality of radiation images of a single object respectively with a plurality of kinds of 
radiation having different energy distributions, different images of at least part of the object being 
embedded in the plurality of the radiation Images, 

ii) obtaining a plurality of image signals, each of which is made up of a series of image signal 
70 components, from the plurality of the radiation images, 

iii) weighting the image signal components of the plurality of the image signals, which Image signal 
components represent corresponding picture elements in the plurality of the radiation images, with 
predetermined weight factors, and 

iv) subtracting the weighted image signal components of the plurality of the image signals from one 
75 another, which image signal components represent corresponding picture elements in the plurality of the 

radiation images, a difference signal representing an image of a specific structure of the object being 
thereby obtained, 

wherein the difference signal is obtained by: 
a) carrying out the convolution of at least a single desired image signal, which is among the plurality of 
20 the image signals, with a mask filter, which has frequency characteristics capable of keeping the signal- 
to-noise ratio of the difference signal high and altering the response characteristics of the difference 
signal with respect to a desired frequency band when the difference signal is obtained by subtracting an 
image signal, that is obtained from the convolution of the desired image signal, and the other image 
signals from one another, a processed image signal being thereby obtained from the convolution, and 
25 b) subtracfing the processed image signal and the other image signals ft'om one another. 

As in the twelfth superposifion processing method for a radiafion image in accordance with the present 
invention, the ninth energy subtraction processing method in accordance with the present invention should 
preferably be modified such that each of the plurality of the image signals may be obtained by scanning a 
sheet-like recording medium, on which the radiation image has been recorded, with a light beam in the 
30 main scanning direction and the sub-scanning direction, the radiation image being thereby read out 
photoelectrically, and 

the processed image signal may be obtained by carrying out the convolution of the desired image 
signal with the mask filter, which has different frequency characteristics with respect to the main scanning 
direction and the sub-scanning direction, such that the frequency response characteristics of the difference 
35 signal with respect to the main scanning direction and the frequency response characteristics of the 
difference signal with respect to the sufc>-scanning direction may become approximately identical with each 
other. 

Also, in the ninth energy subtraction processing method in accordance with the present invention, the 
processing with the convolution may be carried out on all of the plurality of the image signals. 
40 In cases where the processing with the convolution is carried out on all of the plurality of the image 
signals, the term "other image signals" as used herein also means the image signals obtained from the 
convolution. 

With the twelfth superposition processing method for a radiation image in accordance with the present 
invention, the convolution of at least a single desired image signal is canried out with the mask filter, which 

45 has frequency characteristics capable of keeping the signal-to-noise ratio of the addition signal high and 
altering the response characteristics of the addition signal with respect to a desired frequency band. The 
processed image signal obtained from the convolufion and the other image signals are then added to one 
another. Therefore, in the obtained addition signal, noise components can be reduced, and the response 
characteristics with respect to a desired frequency band can be altered. Accordingly, the radiation image 

50 represented by the addition signal can have good image quality and little noise. Further, in the radiation 
image, the desired frequency band has been altered. As a result, a superposition image can be obtained 
which has good image quality and can serve as an effective tool in, particularly, the efficient and accurate 
diagnosis of an illness. In addition, with the twelfth superposition processing method for a radiation image in 
accordance with the present invenfion. the processing with the convolution is carried out on the entire 

55 image signal. Therefore, it is not necessary to carry out a frequency transform, such as the wavelet 
transform or the Fourier transform. Also, it is not necessary to carry out filtering processing on the addition 
signal. Accordingly, the amount of calculation can be kept small, and the apparatus for carrying out the 
superposition processing method for a radiation image in accordance with the present invention can be kept 
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simple. As a result, a superposition image having good image quality can be obtained quickly and at a low 
cost. 

Also, In cases where the processing with the convolution is carried out on all of the plurality of the 
image signals, an addition signal having better image quality can be obtained. 

5 Further, as described above, each of the plurality of the image signals may be obtained by scanning 
the sheet-like recording medium, on which the radiation image has been recorded, with a light beam in two- 
dimensional directions. The processed Image signal may then be obtained by carrying out the convolution 
of the desired image signal with the mask filter, which has different frequency characteristics with respect to 
the main scanning direction and the sub-scanning direction, such that the frequency response characteris- 

70 tics of the addition signal with respect to the main scanning direction and the frequency response 
characteristics of the addition signal with respect to the sub-scanning direction may become approximately 
identical with each other. In such cases, in the radiation image represented by the addition signal, the 
frequency characteristics with respect to the main scanning direction and the frequency characteristics with 
respect to the sub-scanning direction become well-balanced. 

75 Accordingly; from the addition signal, a reproduced image can be obtained which has better image 
quality and can serve as a more effective too! in, particulariy, the efficient and accurate diagnosis of an 
illness. 

The processing carried out in the twelfth superposition processing method for a radiation image in 
accordance with the present invention can also be applied to the energy subtraction processing in the ninth 
20 energy subtraction processing method in accordance with the present invention. Therefore, with the ninth 
energy subtraction processing method in accordance with the present invention, a difference signal 
representing an image having good image quality and containing little noise can be obtained from the 
subtraction processing, and the response characteristics with respect to a desired frequency band can be 
altered. Accordingly, the radiation image represented by the difference signal can have good image quality 
25 and little noise. Further, in the radiation image, the desired frequency band has been altered. As a result, a 
subtraction image can be obtained which has good image quality and can serve as an effective tool in, 
particulariy, the efficient and accurate diagnosis of an illness. In addition, with the ninth energy subtraction 
processing method in accordance with the present invention, the processing with the convolution is carried 
out on the entire image signal. Therefore, it is not necessary to carry out a frequency transform, such as the 
30 wavelet transform or the Fourier transform. Also, it is not necessary to carry out filtering processing on the 
difference signal. Accordingly, the amount of calculation can be kept small, and the apparatus for carrying 
out the energy subtraction processing method in accordance with the present invention can be kept simple. 
As a result, a subtraction image having good image quality can be obtained quickly and at a low cost. 

The present invention further provides a thirteenth superposition processing method for a radiation 
35 image, comprising the steps of: 

i) obtaining a plurality of image signals, which represent a radiation image of a single object or radiation 
images of the single object having been fonmed by delivering radiation to the single object, and which 
have different frequency characteristics, each of the image signals being made up of a series of image 
signal components, and 

40 ii) adding the image signal components of the plurality of the image signals to one another, which image 
signal components represent corresponding picture elements, an addition signal being thereby obtained, 
wherein the addition signal is obtained by: 

a) calculating the dose of radiation delivered to the object, 

b) setting a mask filter for at least a single desired image signal, which is among the plurality of the 
45 image signals, in accordance with the calculated dose of radiation, the mask filter having frequency 

characteristics capable of keeping the signal-to-noise ratio of the addition signal high when the addition 
signal is obtained by adding an image signal, that is obtained from the convolution of the desired Image 
signal, and the other image signals to one another, 

c) carrying out the convolution of the desired image signal with the mask filter, and 

50 d) adding the image signal, which has been obtained from the convolution, and the other image signals 
to one another. 

In the thirteenth superposition processing method for a radiation image in accordance with the present 
invention, the processing with the convolution may be carried out on all of the plurality of the image signals. 
In such cases, the sum of the frequency characteristics of the mask filters employed for the plurality of the 
55 image signals may be equal to 1 at an arbitrary frequency. 

In cases where the processing with the convolution is carried out on all of the plurality of the image 
signals, the term "other image signals" as used herein also means the image signals obtained from the 
convolution. 
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The thirteenth superposition processing method for a radiation image in accordance with the present 
Invention may be modified such that the dose of radiation may be calculated for each of portions of the 
object, the patterns of which are embedded in the radiation Image. 

the mask filter may be set for each of the portions of the object in accordance with the calculated dose 
5 of radiation, and 

the convolution of the desired image signal may be carried out with the mask filter, which has been set 
for each of the portions of the object. 

Further, the Image processing with the convolution may be carried out on each of the Image signals by 
using a single mask filter. 

10 The present invention still further provides a tenth energy subtraction processing method comprising 
the steps of: 

i) fonming a plurality of radiation images of a single object respectively with a plurality of kinds of 
radiation having different energy distributions, different images of at least part of the object being 
embedded in the plurality of the radiation images, 
75 ii) obtaining a plurality of image signals, each of which is made up of a series of image signal 
components, from the plurality of the radiation images. 

iii) weighting the image signal components of the plurality of the Image signals, which image signal 
components represent corresponding picture elements in the plurality of the radiation images, with 
predetemnined weight factors, and 

20 iv) subtracting the weighted Image signal components of the plurality of the image signals from one 
another, which image signal components represent corresponding picture elements in the plurality of the 
radiation images, a difference signal representing an image of a specific structure of the object being 
thereby obtained, 

wherein the difference signal is obtained by: 

25 a) calculating the dose of radiation delivered to the object, 

b) setting a mask filter for at least a single desired Image signal, which is among the plurality of the 
image signals, in accordance with the calculated dose of radiation, the mask filter having frequency 
characteristics capable of keeping the signal-to-noise ratio of the difference signal high when the 
difference signal is obtained by subtracting an Image signal, that is obtained from the convolution of the 

30 desired image signal, and the other image signals from one another, 

c) carrying out the convolution of the desired Image signal with the mask filter, and 

d) subtracting the image signal, which has been obtained from the convolution, and the other image 
signals from one another. 

As In the thirteenth superposition processing method for a radiation image in accordance with the 
35 present invention, the tenth energy subtraction processing method in accordance with the present invention 
may be modified such that the dose of radiation may be calculated for each of portions of the object, the 
patterns of which are embedded in the radiation image, 

the mask filter may be set for each of the portions of the object In accordance with the calculated dose 
of radiation, and 

40 the convolution of the desired Image signal may be carried out with the mask filter, which has been set 
for each of the portions of the object. 

In the tenth energy subtraction processing method in accordance with the present invention, the 
processing with the convolution may be carried out on all of the plurality of the image signals. In such 
cases, the sum of the frequency characteristics of the mask filters employed for the plurality of the Image 
45 signals may be equal to 1 at an arbitrary frequency. 

In cases where the processing with the convolution is carried out on all of the plurality of the image 
signals, the term "other image signals" as used herein also means the image signals obtained from the 
convolution. 

With the thirteenth superposition processing method for a radiation image In accordance with the 
50 present invention, the dose of radiation delivered to the object is calculated. The mask filter to be used for 
the convolution of at least a single desired Image signal is set in accordance with the calculated dose of 
radiation. Therefore, the mask filter, which can yield the addition signal representing a radiation image 
having good image quality, can be set in accordance with the dose of radiation delivered to the object. As a 
result, a superposition image having good Image quality can be obtained regardless of the dose of radiation 
55 delivered to the object. In addition, with the thirteenth superposition processing method for a radiation 
Image In accordance with the present invention, the processing for changing the frequency characteristics is 
carried out on the entire Image signal. Therefore, it is not necessary to canry out a frequency transform, 
such as the wavelet transform or the Fourier transform. Accordingly, the amount of calculation can be kept 
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sma!!, and the apparatus for carrying out the superposition processing method for a radiation innage in 
accordance with the present invention can be kept simple. As a result, a superposition image having good 
Image quality can be obtained quickly and at a low cost. 

Also, in cases where the processing with the convolution is carried out on ail of the plurality of the 
5 image signals, an addition signal having better image quality can be obtained. 

Further, in cases where the sum of the frequency characteristics of the mask filters, which are 
employed for the plurality of the image signals, is equal to 1 at an arbitrary frequency, when the image 
signals obtained from the processing with the convolution are added to one another, it becomes unnec- 
essary for the weighting process to be carried out such that the addition ratio of the image signals may be 
^ 10 equal to 1. Therefore, the operation time can be kept short, and the addition process can be carried out 
quickly. 

Furthermore, as described above, the dose of radiation may be calculated for each of portions of the 
object, the pattems of which are embedded in the radiation image. The mask filter may then be set for each 
of the portions of the object in accordance with the calculated dose of radiation. In such cases, it becomes 

75 possible to obtain an addition signal representing a radiation image which has good image quality with 
respect to each of the portions of the object and can serve as an effective tool in, particularly, the efficient 
and accurate diagnosis of an Illness at each of the portions of the object. In this manner, a superposition 
image having better image quality can be obtained. 

Moreover, in cases where the image processing with the convolution is carried out on each of the 

20 image signals by using a single mask filter, the number of the mask filters, which are to be stored in an 
apparatus for carrying out the superposition processing method for a radiation image in accordance with the 
present invention, can be kept small. Therefore, the apparatus can be kept simple. 

The processing carried out in the thirteenth superposition processing method for a radiation image in 
accordance with the present invention can also be applied to the energy subtraction processing in the tenth 

25 energy subtraction processing method in accordance with the present invention. Therefore, with the tenth 
energy subtraction processing method in accordance with the present invention, a difference signal 
representing an image having good image quality and containing little noise in accordance with the dose of 
radiation delivered to the object can be obtained from the subtraction processing. Also, the image 
processing, which changes the frequency characteristics of the image signal, is carried out on the entire 

30 image signal. Therefore, it is not necessary to carry out a frequency transform, such as the wavelet 
transform or the Fourier transform. Accordingly, the amount of calculation can be kept small, and the 
apparatus for carrying out the tenth energy subtraction processing method in accordance with the present 
invention can be kept simple. As a result, a subtraction image having good image quality can be obtained 
quickly and at a low cost. 

35 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a schematic view showing how radiation images are recorded on stimulable phosphor sheets 
in an embodiment of the superposition processing method for a radiation image in accordance with the 
40 present invention. 

Figure 2 is a perspective view showing an apparatus for reading out a radiation image from a stimulable 
phosphor sheet, on which the radiation image has been stored, 

Figure 3 is a block diagram showing an apparatus for carrying out a first embodiment of the 
superposition processing method for a radiation image in accordance with the present invention, 

45 Figure 4 is a flow chart showing how a wavelet transform is carried out. 

Figure 5 is an explanatory view showing images represented by wavelet transform factor signals, 
Figures 6A and 6B are graphs showing modulation transfer functions (MTF's) of image signals. 
Figures 7A and 78 are graphs showing Winer spectra of image signals. 
Figures 8A, SB, and 8C are graphs showing DQE's for different frequency bands, 

50 Figure 9 is a graph showing a weight table, 

Figure 10 is a flow chart showing how an inverse wavelet transform is carried out, 

Figure 1 1 is a flow chart showing how a sub-band transform is carried out. 

Figure 12 is an explanatory view showing filters for carrying out a Fourier transform. 

Figure 13 is an explanatory view showing images, which are represented by transform factor signals 

55 obtained from a wavelet transform or a sub-band transform, 

Figure 14 is a schematic view showing how a radiation image is recorded on a single stimulable 
phosphor sheet. 
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Figure 15 is a perspective view showing an apparatus for detecting light ennitted by two surfaces of a 
stimulable phosphor sheet, 

Figure 16 is a perspective view showing a different apparatus for detecting light emitted by two surfaces 
of a stimulable phosphor sheet, 

5 Figure 1 7 is a block diagram showing how frequency processing is carried out on analog output signals, 
Figure 18 is a graph showing a filter for carrying out frequency processing on an output signal SA, 
Figure 19 is a graph showing a filter for carrying out frequency processing on an output signal SB, 
Figure 20 is a schematic view showing how radiation images, which are to be subjected to energy 
subtraction processing, are recorded on stimulable phosphor sheets, 

10 Figure 21 Is a block diagram showing an apparatus for carrying out a first embodiment of the energy 
subtraction processing method in accordance with the present invention, 
Figure 22 is a graph showing a basic wavelet function employed in a wavelet transform. 
Figure 23 is a diagram showing signals, which are obtained by carrying out a wavelet transform on an 
original signal Sorg and then carrying out an inverse wavelet transform for each of frequency bands, 

15 Figure 24 is a diagram showing signals, which are obtained by canying out a Fourier transform on the 
original signal Sorg and then carrying out an inverse Fourier transform for each of the frequency bands. 
Figure 25 Is a block diagram showing an apparatus for carrying out a second embodiment of the 
superposition processing method for a radiation image in accordance with the present invention, 
Figure 26 is an explanatory view showing how an unsharp mask signal is generated, 

20 Figure 27 is a graph showing frequency characteristics of a filter used for an image signal S2, 
Figure 28 is a graph showing frequency characteristics of a filter used for an image signal SI, 
Figure 29 is a graph showing frequency characteristics of a filter F3 used for an image signal 81, 
Figure 30 is a graph showing frequency characteristics of a filter F4 used for an image signal S2, 
Figure 31 is a graph showing frequency characteristics of the sum of the filters F3 and F4, 

25 Figure 32 is a block diagram showing an apparatus for carrying out a second embodiment of the energy 
subtraction processing method in accordance with the present invention, 

Figure 33 is a block diagram showing an apparatus for carrying out a third embodiment of the 
superposition processing method for a radiation image in accordance with the present invention, 
Figure 34 is a graph showing frequency characteristics of a filter F2 used for the convolution of an image 
30 signal S2, 

Figure 35 is a graph showing frequency characteristics of a filter F3 used for the convolution of an 
addition signal, 

Figure 36 is a graph showing frequency characteristics of a filter FV used for the convolution of an 
image signal Si, 

35 Figure 37 is a graph showing frequency characteristics of a filter F2* used for the convolution of the 
image signal S2, 

Figure 38 is a graph showing frequency characteristics of the sum of the filters F1' and F2\ 
Figure 39 is a block diagram showing an apparatus for carrying out a third embodiment of the energy 
subtraction processing method in accordance with the present invention, 
40 Figure 40 is a perspective view showing an apparatus for reading out a radiation image from a stimulable 
phosphor sheet, on which the radiation image has been stored, 

Figure 41 is a block diagram showing an apparatus for carrying out a fourth embodiment of the 
superposition processing method for a radiation image in accordance with the present invention, 
Figure 42 is a graph showing frequency characteristics of a filter used for the convolution of an image 
45 signal S2 (radiation dose: 0,1 mR), 

Figure 43 is a graph showing frequency characteristics of a filter used for the convolution of an image 
signal S2 (radiation dose: I.OmR). 

Figure 44 is a graph showing frequency characteristics of a filter used for the convolution of an image 
signal S2 (radiation dose: lOmR), and 
50 Figure 45 is a block diagram showing an apparatus for carrying out a fourth embodiment of the energy 
subtraction processing method in accordance with the present invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

55 The present invention will hereinbelow be described in further detail with reference to the accompany- 
ing drawings. 

Rgure 1 shows how radiation 2, which has passed through a single object 1, is inradiated to two 
stimulable phosphor sheets 4A and 4B, 
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As illustrated in Figure 1 , the first stimulable phosphor sheet 4A and the second stimulable phosphor 
sheet 4B are superposed one upon the other, and a radiation source 3 is activated to produce the radiation 
2. The radiation 2, which has been produced by the radiation source 3, passes through the object 1 . The 
radiation 2, which has passed through the object 1, impinges upon the first stimulable phosphor sheet 4A 

5 and the second stimulable phosphor sheet 4B. In this manner, radiation images of the object 1 are stored 
on the first stimulable phosphor sheet 4A and the second stimulable phosphor sheet 4B. 

Thereafter, the radiation images are read out from the first stimulable phosphor sheet 4A and the 
second stimulable phosphor sheet 48 by using an Image read-out means shown in Figure 2, and image 
signals representing the radiation images are thereby obtained. Specifically, the first stimulable phosphor 

70 sheet 4A is moved by a sub-scanning means 9 in the sub-scanning direction indicated by the arrow Y. The 
sub-scanning means 9 may be constituted of an endless belt, or the like. At the same time, a laser beam 
11, which serves as stimulating rays, is produced by a laser beam source 10. The laser beam 11 is 
deflected by a scanning mirror 12 and caused to scan the stimulable phosphor sheet 4A in the main 
scanning directions indicated by the double-headed arrow X. When the stimulable phosphor sheet 4A is 

75 exposed to the laser beam 1 1 , it emits light 1 3 in proportion to the amount of energy stored thereon during 
its exposure to the radiation 2. The emitted light 13 enters a light guide member 14, which is made fi'om a 
transparent acrylic plate, from its one edge face. The emitted light 13 is guided through repeated total 
reflection inside of the light guide member 14 and detected by a photomultiplier 15. The photomultiplier 15 
generates an output signal SA corresponding to the amount of the emitted light 13, i.e. representing the 

20 radiation image stored on the stimulable phosphor sheet 4A. 

The output signal SA is logarithmically amplified by a logarithmic amplifier 16 and is then converted by 
an analog-to-digital converter 17 into a digital image signal SI. The digital image signal Si is stored on a 
storage medium 18, such as a magnetic disk. Thereafter, the radiation image stored on the second 
stimulable phosphor sheet 48 is read out in the same manner as that described above, and an output signal 

25 SB representing the radiation image is thereby obtained. The output signal SB is logarithmically amplified 
by the logarithmic amplifier 16 and is then converted by the analog-to-digitat converter 17 into a digital 
image signal S2. The digital image signal S2 is stored on the storage medium 18. 

Thereafter, a superposition processing is carried out on the image signals SI and S2. Figure 3 is a 
block diagram showing an apparatus for carrying out a first embodiment of the superposition processing 

30 method for a radiation image in accordance with the present invention. First, the image signals SI and S2 
are read from an image file 18A and an image file 18B in the storage medium 18 and fed into a wavelet 
transform means 19, The wavelet transform means 19 carries out a wavelet transform of each of the two 
image signals SI and S2 and decomposes each of the two image signals S1 and S2 into a plurality of 
wavelet transfomn factor signals, each of which is of one of a plurality of different fi'equency bands. How the 

35 wavelet transform is carried out will be described hereinbelow. 

Figure 4 is a flow chart showing how the wavelet transform is carried out on each of the image signals 
SI and S2. As an aid in facilitating the explanation, how the wavelet transform of the image signal SI is 
carried out will be described hereinbelow. 

In this embodiment, an orthogonal wavelet transform, in which the respective wavelet transform factors 

40 are orthogonal, is carried out. The orthogonal wavelet transform is described in the aforesaid literature of 
Marc Antonini, et al. 

As illustrated in Rgure 4, filtering processing is carried out with a function g and a function h, which are 
obtained from the basic wavelet function, on the image signal components of the image signal S1 
representing picture elements in the radiation image, which are located along the main scanning direction. 

45 Specifically, the filtering processing on the image signal components of the image signal Si representing 
each row of the picture elements, which are arrayed along the main scanning direction, is carried out with 
the function g and the function h each time the position of the filtering processing is shifted by a single 
picture element in the sub-scanning direction. In this manner, wavelet transform factor signals WgO and 
Who with respect to the main scanning direction of the image signal SI are obtained. 

50 The function g and the fijnction h can be uniquely obtained from the basic wavelet function. For 
example, the function h has tiie characteristics shown in Table 1 below. In Table 1 , a function h' is the one 
which is used when an inverse wavelet transform is carried out on an image signal having been subjected 
to the wavelet transform. As will be understood from Formula (4), the function g can be obtained from the 
function h'. and a function g* to be used during the inverse wavelet transform can be obtained from the 

55 function h. 
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The wavelet transform factor signals WgO and WhO are obtained in the manner described above. 
Thereafter, for each of the wavelet transform factor signals WgO and WhO, the signal components 
representing the picture elements located along the main scanning direction are thinned out alternately. In 
this manner, the number of the picture elements located along the main scanning direction is reduced to 
1/2. Filtering processing is then can'ied out with the function g and the function h on the signal components 
of the thinned-out wavelet transform factor signals WgO and WhO, which signal components represent 
picture elements located along the sub-scanning direction. From the filtering processing, wavelet transform 
factor signals WWuo. WVyo, VWuo, and Wyo are obtained. 

Thereafter, for each of the wavelet transform factor signals WWuo, WVuo. VWuo, and Wuo. the signal 
components representing the picture elements located along the sub-scanning direction are thinned out 
alternately. In this manner, the number of the picture elements located along the sub-scanning direction is 
reduced to 1/2. As a result, the number of the picture elements represented by each of the wavelet 
transform factor signals WWuo, WVuo, VWuo, and Wuo becomes equal to 1/4 of the number of the picture 
elements represented by the image signal SI . Rltering processing is then carried out with the function g 
and the function h on the signal components of the wavelet transform factor signal Wuo, which represent 
picture elements located along the main scanning direction. 

Specifically, the filtering processing on the image signal components of the wavelet transform factor 
signal Wuo representing each row of the picture elements, which are an'ayed along the main scanning 
direction, is carried out with the function g and the function h each time the position of the filtering 
processing Is shifted by a single picture element In the sub-scanning direction. In this manner, wavelet 
transform factor signals Wgl and Whi with respect to the main scanning direction of the wavelet transfonm 
factor signal Wuo are obtained. 

The number of the picture elements represented by the wavelet transform factor signal Wuo is equal to 
1/2 of the number of the picture elements, which are represented by the original image signal, both In the 
main scanning direction and in the sub-scanning direction. Therefore, the resolution of the image repre- 
sented by the wavelet transform factor signal Wuo 's equal to 1/2 of the resolution of the image represented 
by the original image signal. Accordingly, as a result of the filtering processing carried out with the function 
g and the function h on the wavelet transform factor signal Wuo, the wavelet transform factor signals Wgl 
and Whi representing the frequency components, which are lower than the frequency components 
represented by the wavelet transfomn factor signal Wuo and which are among the frequency components of 
the original image signal, are obtained. 

The wavelet transform factor signals Wgl and Whi are obtained in the manner described above. 
Thereafter, for each of the wavelet transfonm factor signals Wgl and Whi, the signal components 
representing the picture elements located along the main scanning direction are thinned out alternately. In 
this manner, the number of the picture elements located along the main scanning direction is reduced even 
further to 1/2. Filtering processing is then canied out with the function g and the function h on the signal 
components of the thinned-out wavelet transform factor signals Wgl and Whi, which signal components 
represent picture elements located along the sub-scanning direction. From the filtering processing, wavelet 
transform factor signals WWyi, WVui, WVui, and Wui are obtained. 

Thereafter, for each of the wavelet transform factor signals WWyi, WVui, WVui, and Wui, the signal 
components representing the picture elements located along the sub-scanning direction are thinned out 
alternately. In this manner, the number of the picture elements located along the sub-scanning direction is 
reduced to 1/2. As a result, the number of the picture elements represented by each of the wavelet 
transform factor signals WWui, WVui, WVui, and Wui becomes equal to 1/16 of the number of the picture 
elements represented by the image signal Si . 

Thereafter, in the same manner as that described above, filtering processing is carried out with the 
function g and the function h on the signal components of the thinned-out wavelet transform factor signal 
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Wui. which represent picture elements located along the main scanning direction. Then, for each of the 
wavelet transform factor signals, which have thus been obtained, the signal components representing the 
picture elements located along the main scanning direction are thinned out alternately. Rltering processing 
is then carried out with the function g and the function h on the signal components of the thinned-out 
5 wavelet transform factor signals, which signal components represent picture elements located along the 
sub-scanning direction. From the filtering processing, wavelet transform factor signals WWu2. WVua. VWu2. 
and VVu2 are obtained. 

The wavelet transform described above is iterated N number of times, and wavelet transform factor 
signals \A/Wuo through WW^n. WVuo through WVun, VWuq through VWun. and Wun are thereby obtained. 

70 The number of the picture elements represented by each of the wavelet transform factor signals WWun, 
WVgN. VWuN. and Wun. which are obtained from the N'th wavelet transfonm, is equal to (1/2)" of the 
number of the picture elements, which are represented by the original image signal, both in the main 
scanning direction and in the sub-scanning direction. Therefore, as the value of N becomes larger, each 
wavelet transform factor signal Is of a lower frequency band and represents lower frequency components 

75 among the frequency components of the original image signal. 

Accordingly, a wavelet transform factor signal WWui (wherein 1=0 to N, this also applies to the 
descriptions below) represents a change in the frequency of the image signal SI both in the main scanning 
direction and in the sub-scanning direction. As the value of i becomes larger, the wavelet transform factor 
signal WWui becomes a lower frequency signal. A wavelet transform factor signal WVui represents a 

20 change in the frequency of the image signal S1 in the main scanning direction. As the value of i becomes 
larger, the wavelet transform factor signal WVui becomes a lower frequency signal. A wavelet transform 
factor signal NA/Vji represents a change in the frequency of the image signal Si in the sub-scanning 
direction. As the value of i becomes larger, the wavelet transform factor signal VWyi becomes a lower 
frequency signal. 

25 Figure 5 shows images represented by the wavelet transform factor signals for the respective frequency 
bands. As an aid in facilitating the explanation, only the images represented by the wavelet transform factor 
signals obtained from the tirst, second, and third wavelet transfomns are shown in Figure 5. In Figure 5, the 
wavelet transform factor signal WWus represents an image, which is obtained by reducing the original 
image to (1/2)^ in each of the main scanning direction and the sub-scanning direction. 

30 The wavelet transforms are carried out on the image signal S2 in the same manner as that described 
above. From the wavelet transforms, wavelet transform factor signals WWuo through WWln. WVlo through 
WVtN, VWlo through VWuj, and Wln are thereby obtained for the respective frequency bands. 

The wavelet transform factor signals, which have thus been obtained by carrying out the wavelet 
transforms on the image signals Si and S2, are fed into a weighting and superposition means 20. In the 

35 weighting and superposition means 20, a weighting process is carried out such that tiie value of the weight 
factor with respect to the frequency band, which has a low signal-to-noise ratio, may be rendered smaller 
tiian the value of the weight factor with respect to the frequency band, which has a high signal-to-noise 
ratio. How the value of the weight factor is determined will be described below. 

The image signals Si and S2, which have been obtained from tiie two stimulable phosphor sheets 4A 

40 and 4B, respectively have the modulation transfer functions (MTF's, i.e. frequency dependency characteris- 
tics) shown in Figures 6A and 6B. The MTF can be obtained by recording a contrast transfer function chart 
(CTF chart) and represents the level of the resolution of the image signal with respect to each frequency 
band. Specifically, as illustrated in Figure 6A, the MTF1 of tiie image signal S1, which has been obtained 
from the stimulable phosphor sheet 4A located at the position closer to the radiation source during the 

45 image recording operation, takes a large value up to the high frequency band. Thus, the image signal S1 
carries the information up to the high frequency band. On the other hand, as illustrated in Figure 6B, the 
MTF2 of the image signal S2, which has been obtained from the stimulable phosphor sheet 48 located at 
the position remote from the radiation source during the image recording operation, takes a smaller value 
on the side of the high frequency band than the MTF1 of the image signal SI. Thus, in the image signal S2, 

50 the amount of information in the high frequency band is small. This indicates that the information in the high 
frequency band of the image signal S2 contains noise, due to scattered radiation during the image 
recording operation, or the like, and that the fine information on the high frequency band side has been 
rendered unsharp due to the location of the stimulable phosphor sheet 48 remote from the radiation source. 
Therefore, the wavelet transform factor signals, which are of a single frequency band, are weighted with the 

55 weight factors, such that the values of the weight factors may be varied in accordance with the MTF. The 
weighted wavelet transform factor signals, which are of a single frequency band, are then added to one 
another. How the values of the weight factors for the wavelet transform factor signals are determined will be 
described herelnbelow. 
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The frequency characteristics MTFl and MTF2 of the image signals shown in Figures 6A and 6B are 
calculated. Also, as illustrated in Figures 7A and 7B, frequency characteristics Winer 1 and Winer 2 of noise 
of the image signals are calculated. Each of the Winer 1 and the Winer 2 represents the amount of 
dispersion of the noise image signal, which has been obtained by recording an image of only the noise, i.e. 
by carrying out the image recording operation without the object lying, with respect to each frequency. 
Specifically, as for the Winer 1, the image of only the noise is recorded, and a noise image signal Image 
(X1) is obtained from the upper stimulable phosphor sheet 4A. The calculation with Formula (5) is carried 
out on the Image (X1) 



RMS2 = /{image {X)}2 dx-image pQZ*'"^^"* (5) 

and the value ot RMS^ is thus obtained. The value of RMS^ is plotted for each frequency. In this manner, 
the Winer 1 shown in Figure 7A is obtained. In the same manner, the Winer 2 shown in Figure 7B is 
obtained. 

A DQE index is defined by Fonmula (6). 

DQE oc (MTF)2/Winer (6) 

Formula (6) indicates that a higher DQE value represents better image quality. The DQE is calculated for 
each frequency. 

Thereafter, the image signals, Image 1 (X) and Image 2 (X), for each frequency band, which are 
obtained when the MTFl and MTF2 are obtained, are added to each other, and an addition image signal 
add(t) is thereby obtained. The calculation is carried out with Formula (7). 

add{t) = t X Image 1 (X) + (1-t) x Image 2 (X) (7) 

The value of t is changed between 0 and 1, and a plurality of addition image signals add(t) are thereby 
obtained. The DQE is calculated for each of the addition image signals add(t) and plotted on the graph, in 
which the value of t is plotted on the horizontal axis, and the DQE is plotted on the vertical axis. Figures 8A, 
88, and 8C are graphs showing the relationship between t and DQE having been obtained for each of the 
plurality of frequency bands. As illustrated in Rgure 8A. when the frequency band is 1 cycle/mm (indicated 
as Ic/mm in Figure 8A), the DQE takes the largest value at t = 0.5. Also, as illustrated in Figure 88, when 
the frequency band is 2 cycles/mm, the DQE takes the largest value at t = 0.7. Further, as illustrated in 
Rgure 8C, when the frequency band is 3 cycles/mm, the DQE takes the largest value at t = 0.9. 

When the value of t, which is associated with the largest value of DQE with respect to each of the 
frequency bands, is plotted, the weight table shown in Figure 9 can be obtained. The wavelet transform 
factor signals, which are of a single frequency band, are weighted in accordance with the weight table 
shown in Figure 9, and the weighted wavelet transform factor signals, which are of a single frequency band, 
are then added to one another. 

Specifically, the weighted additions of the wavelet transform factor signals, which are of a single 
frequency band, are carried out with Formula (8) 

WW) = t • WWui + (1 - 1) WWli (8) 
WV| = t • WVui + (1 • t) WVu 
VW, = t • VWui + (1 - 1) VWu 
W, = t • Wui + (1 - 1) Wli 

For example, the wavelet transform factor signal WW^ contains more noise and a smaller amount of 
information than the wavelet transform factor signal WWui. Therefore, when an addition wavelet transform 
factor signal WWi for the high frequency band is to be obtained, the value of t is set to be large. 
Specifically, the addition wavelet transform factor signal WWi is calculated with Formula (9). 

WWi = 0.8 X WWui + 0.2 X WWli (9) 

When addition wavelet transform factor signals WVi and VWi are to be obtained, the weighting may be 
carried out in the same manner as that for the addition wavelet transfonm factor signal WWi . 

Also, the differences in the amount of noise and the amount of infonmation between the wavelet 
transform factor signals WWl2 and WWu2 are no so large as the differences between the wavelet transform 
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factor signals WWli and WWyi. Therefore, in cases where an addition wavelet transform factor signal WW2 
for a frequency band lower than the frequency band of the addition wavelet transform factor signal WWi is 
to be obtained, the value of t is set to be approximately 0.6. The addition wavelet transfomn factor signal 
WW2 is thus calculated with Fomnula (10). 

5 

WW2 = 0.6 X WWu2 + 0.4 X WWl2 (10) 

When addition wavelet transform factor signals WV2 and VW2 are to be obtained, the weighting may be 
carried out in the same manner as that for the addition wavelet transform factor signal WW2. 
70 Further, the wavelet transform factor signals WWl3 and WWu3 represent approximately the same 
amounts of information. Therefore, in cases where an addition wavelet transform factor signal WW3 for a 
frequency band lower than the frequency band of the addition wavelet transform factor signal WW2 is to be 
obtained, the value of t is set to be 0.5. The addition wavelet transform factor signal WW3 is thus calculated 
with Formula (11). 

75 

WW3 = 0.5 X WWu3 + 0.5 X WWl3 (11) 

When addition wavelet transform factor signals WVa and VW3 are to be obtained, the weighting may be 
carried out in the same manner as that for the addition wavelet transform factor signal WW3. 

20 Furthemnore, the wavelet transform factor signals WWl4 and WWu4, the wavelet transform factor signals 

WWl5 and WWus the wavelet transform factor signals WWln and WWun respectively represent 

approximately the same amounts of information. Therefore, in cases where addition wavelet transform factor 
signals WW4. WW5, WWn for frequency bands lower than the frequency band of the addition wavelet 
transform factor signal WW3 are to be obtained, the value of t Is set to be 0.5. 

25 The values of the weight factors are determined in the manner described above. Therefore, regardless 
of the characteristics of the MTF and the Winer of the original image, appropriate values of the weight 
factors can be determined for each frequency. 

In the manner described above, the addition wavelet transform factor signals WWi through WWn, WVi 
through WVn, VWi through VWn, and Wi through Wr, are obtained in the weighting and superposition 

30 means 20. Thereafter, In an inverse wavelet transform means 21, an inverse wavelet transform is carried out 
on each of the addition wavelet transform factor signals. How the inverse wavelet transform is carried out 
will be described hereinbelow. 

Figure 10 is a flow chart showing how an inverse wavelet transform is carried out. 

As illustrated in Figure 10, each of the addition wavelet transfomn factor signals . Wn, VWn. WVn, and 

35 WWn is subjected to the processing for leaving a space, which has a length equal to the length of a single 
picture element, between adjacent picture elements located along the sub-scanning direction. (In Figure 10. 
this processing is expressed as x2.) Filtering processing Is then earned out with a function h\ which is 
different from the aforesaid function h, on the signal components of the addition wavelet transform factor 
signal Wn provided with the spaces, which signal components represent picture elements located along the 

40 sub-scanning direction. Also, filtering processing is carried out with a function g*. which is different from the 
aforesaid function g, on the signal components of the addition wavelet transform factor signal VWn provided 
with the spaces, which signal components represent picture elements located along the sub-scanning 
direction. 

Specifically, the filtering processing on the image signal components of the addition wavelet transform 
45 factor signal Wn representing each column of the picture elements, which are arrayed along the sub- 
scanning direction, is carried out with the function h' each time the position of the filtering processing is 
shifted by a single picture element in the main scanning direction. Also, the filtering processing on the 
image signal components of the addition wavelet transform factor signal VWn representing each column of 
the picture elements, which are arrayed alorig the sub-scanning direction, is canned out with the function g' 
50 each time the position of the filtering processing is shifted by a single picture element in the main scanning 
direction. In this manner, inverse wavelet transform factor signals are obtained from the addition wavelet 
transform factor signals Wn and VWn. The inverse wavelet transform factor signals are then doubled and 
added to each other. In this manner, an inverse wavelet transform factor signal WhN' is obtained. 

As described above, the function for the wavelet transform and the function for the Inverse wavelet 
55 transform are different from each other. Specifically, it Is difficult to design functions, which become 
identical in the wavelet transfonm and the inverse wavelet transfomn, i.e. which are the orthogonal functions. 
Therefore, it Is necessary to relieve the conditions of orthogonality, continuity, shortness of function, or 
symmetry. Accordingly, In this embodiment, the conditions of orthogonality are relieved, and the functions 
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satisfying the other conditions are thereby selected. 

As described above. In this embodiment, the functions h and g for the wavelet transform and the 
functions and g* for the Inverse wavelet transform are biorthogonal different functions. Therefore, the 
addition signal of the image signals Si and S2 can be perfectly restored by subjecting the addition wavelet 
5 transform factor signals W], VWi. WVi, and WW| to the inverse wavelet transform with the functions h* and 

g'- 

Also, filtering processing is carried out with the function h' on the signal components of the addition 
wavelet transform factor signal WVn, which represent picture elements located along the sub-scanning 
direction. Also, filtering processing is carried out with the function g' on the signal components of the 
10 addition wavelet transform factor signal WWn, which represent picture elements located along the sub- 
scanning direction. In this manner, Inverse wavelet transform factor signals are obtained from the addition 
wavelet transform factor signals WVn and WWn. The inverse wavelet transform factor signals are then 
doubled and added to each other. In this manner, an Inverse wavelet transform factor signal WgN' is 
obtained. 

75 Thereafter, each of the inverse wavelet transfom factor signals WhN' and WgN' is subjected to the 
processing for leaving a space, which has a length equal to the length of a single picture element, between 
adjacent picture elements located along the main scanning direction. Filtering processing is then carried out 
with the function h* on the signal components of the inverse wavelet transform factor signal WhN\ which 
represent picture elements located along the main scanning direction. Also, filtering processing is carried 

20 out with the function g' on the signal components of the Inverse wavelet transform factor signal WgN', which 
represent picture elements located along the main scanning direction. In this manner, inverse wavelet 
transform factor signals are obtained from the inverse wavelet transform factor signals WhN' and WgN'. The 
inverse wavelet transform factor signals, which have thus been obtained, are then doubled and added to 
each other. In this manner, an addition inverse wavelet transform factor signal Wn-i' is obtained. 

25 Thereafter, each of the addition inverse wavelet transform factor signal Wn-i' and the addition wavelet 
transform factor signals VWn-i, WVn-i. and WWn-i Is subjected to the processing for leaving a space, 
which has a length equal to the length of a single picture element, between adjacent picture elements 
located along the sub-scanning direction. Filtering processing is then carried out with the function h' on the 
signal components of the addition inverse wavelet transform factor signal VVn-i\ which represent picture 

30 elements located along the sub-scanning direction. Also, filtering processing is canied out with the function 
g' on the signal components of the addition wavelet transform factor signal VWn-i. which represent picture 
elements located along the sub-scanning direction. 

Specifically, the filtering processing on the image signal components of the addition inverse wavelet 
transform factor signal Wn-i' representing each column of the picture elements, which are arrayed along 

35 the sub-scanning direction, is carried out with the function h' each time the position of the filtering 
processing Is shifted by a single picture element In the main scanning direction. Also, the filtering 
processing on the image signal components of the addition wavelet transform factor signal VWn-i 
representing each column of the picture elements, which are arrayed along the sub-scanning direction, is 
carried out with the function g' each time the position of the filtering processing is shifted by a single 

40 picture element in the main scanning direction. In this manner, inverse wavelet transform factor signals are 
obtained from the addition inverse wavelet transform factor signal Wn-i' and the addition wavelet transform 
factor signal VWn-i- The inverse wavelet transform factor signals, which have been obtained in this manner, 
are then doubled and added to each other. In this manner, an Inverse wavelet transform factor signal WhN- 
1' is obtained. 

45 Also, filtering processing is carried out with the function h' on the signal components of the addition 
wavelet transform factor signal WVn-i, which represent picture elements located along the sub-scanning 
direction. Also, filtering processing Is carried out with the function g' on the signal components of the 
addition wavelet transform factor signal WWn-i, which represent picture elements located along the sub- 
scanning direction. In this manner, inverse wavelet transform factor signals are obtained from the addition 

50 wavelet transform factor signals WVn-i and WWn-i. The inverse wavelet transform factor signals, which 
have been obtained in this manner, are then doubled and added to each other. In this manner, an inverse 
wavelet transform factor signal WgN-1' is obtained. 

Thereafter, each of the Inverse wavelet transform factor signals WhN-1' and WgN-r is subjected to the 
processing for leaving a space, which has a length equal to the length of a single picture element, between 

55 adjacent picture elements located along the main scanning direction. Filtering processing is then carried out 
with the function h' on the signal components of the Inverse wavelet transform factor signal WhN-1 which 
represent picture elements located along the main scanning direction. Also, filtering processing is carried 
out with the function g' on the signal components of the inverse wavelet transform factor signal WgN-r, 
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which represent picture elements located along the nnain scanning direction. In this manner, inverse wavelet 
transform factor signals are obtained from the Inverse wavelet transform factor signals WhN-V and WgN-r. 
The inverse wavelet transfomn factor signals, which have been obtained in this manner, are then doubled 
and added to each other. In this manner, an addition Inverse wavelet transform factor signal Wn-s' 'S 
5 obtained. 

Thereafter, addition inverse wavelet transform factor signals W/ (wherein i = -1 to N) are sequentially 
created, and an addition inverse wavelet transform factor signal VV-i* is finally obtained. The finally 
obtained addition inverse wavelet transform factor signal W-r serves as an addition image signal Sadd of 
the image signals SI and S2. 

70 The addition inverse wavelet transform factor signal W-i\ which has thus been obtained, is subjected 
to predetermined image processing in an image processing means 22, fed into an image reproducing 
means 23, and used for the reproduction of a visible radiation Image. 

The image reproducing means may be a display device, such as a cathode ray tube (CRT) display 
device, or a recording apparatus for recording an Image on photosensitive film by scanning the photosen- 

75 sitive film with a light beam. 

In the manner described above, the wavelet transform is carried out on each of the two image signals 
SI and S2, and each of the two image signals Si and S2 is thereby decomposed into a plurality of the 
wavelet transform factor signals, each of which is of one of a plurality of different frequency bands. The 
wavelet transform factor signals, which are of a single frequency band, are weighted with the weight factors, 

20 such that the values of the weight factors may be varied for the wavelet transform factor signals of the 
different frequency bands. As for the wavelet transform factor signals for the high frequency band, the value 
of the weight factor with respect to the Image signal, which has been obtained from the stimulable phosphor 
sheet located at the position closer to the radiation source, is set to be larger than the value of the weight 
factor with respect to the image signal, which has been obtained from the stimulable phosphor sheet 

25 located at the position remote from the radiation source. In this manner, it is possible to obtain an addition 
wavelet transfomn factor signal, in which the noise components contained in the image signal obtained from 
the stimulable phosphor sheet located at the position remote from the radiation source have been reduced. 
Each of the addition wavelet transform factor signals, which have thus been obtained, is subjected to the 
inverse wavelet transform, and an addition signal is thereby obtained. From the addition signal thus 

30 obtained, an image having good Image quality and containing little noise component can be reproduced. 

In the embodiment described above, the functions having the characteristics shown in Table 1 are 
employed as the functions h and h' for the wavelet transform. Alternatively, the functions having the 
characteristics shown in Table 2 or Table 3 may be employed as the functions h and h* for the wavelet 
transform. 

35 

Table 2 



40 



45 



50 



n 


0 


i 1 


±2 


t 3 


± 4 


2-1/2^ 
2-1/2^. 


0.6 
17/28 


0.25 
73/280 


-0.05 
-3/56 


0 

-3/280 


0 
0 


Table 3 


n 


0 


i 1 


±2 


- 3 


± 4 


2-1/2^ 
2-i/2h' 


45/64 
1/2 


19/64 
1/4 


-1/8 
0 


-3/64 
0 


3/128 
0 



Also, any of other functions, which can can7 out the wavelet transform, may be employed. For 
example, functions, which are not biorthogonal nor symmetric, but which are orthogonal, may be employed. 

The wavelet transform may be carried out by using the functions, each of which is symmetric with 
55 respect to the axis of n = 0 as shown in Table 1, 2, or 3, or functions, each of which is asymmetric with 
respect to the axis of n = 0. In cases where the wavelet transform has been carried out by using the 
functions, each of which is asymmetric with respect to the axis of n = 0. the inverse wavelet transform is 
carried out by using the functions, v4iich are obtained by inverting the functions used for the wavelet 



26 



EP 0 677 780 A2 

transform with respect to the axis of n = 0. Specifically, the functions g' and h* for carrying out the inverse 
wavelet transform with respect to the functions g and h, each of which is asymmetric with respect to the 
axis of n = 0. are expressed as 

5 g[nl = g*[-n] (12) 

h[nl = h'[-n] 

wherein [-n] represents the inversion with respect to the axis. 

In the embodiment described above, each of the two image signals Si and S2 is decomposed by the 

70 wavelet transform Into a plurality of the wavelet transform factor signals, each of which is of one of a 
plurality of different frequency bands. Alternatively, each of the two image signals S1 and S2 may be 
decomposed by a sub-band transform into a plurality of the transform factor signals, each of which is of one 
of a plurality of different frequency bands. With the wavelet transform, the filtering processing of the image 
signal is carried out with a single kind of the function, and the plurality of the transform factor signals, each 

75 of which is of one of a plurality of different frequency bands, are thereby obtained successively. On the 
other hand, with the sutnband transform, the filtering processing of the image signal is carried out with a 
plurality of functions having different periods, and the*plurality of the transform factor signals, each of which 
is of one of a plurality of different frequency bands, are thereby obtained with a single, simultaneous 
process. 

20 For example, in the embodiment described above, the filtering processing of the image signal is carried 
out with the function g and the function h, and the image signals for the plurality of the frequency bands are 
thereby obtained. On the other hand, with the sub-band transform, as illustrated In Figure 11, the filtering 

processing of the image signal 81 is carried out with a plurality of functions g1, hi, g2, h2 gN. hN such 

that the periods of the function g and the function h become two times, four times, .... 2" times. In this 

25 manner, transform factor signals WWyi. through WWun. WVui through WVyw, \AA/ui through VWun, and 
Wui through Wun are obtained for the plurality of frequency bands. Also, in the same manner, from the 
image signal S2, transfonm factor signals N/VW^ through WWln, WV^ through WVln, VWu through VWln, 
and Wli through Wln are obtained for the plurality of frequency bands. 

Thereafter, in the same manner as that in the aforesaid wavelet transform, the transform factor signals. 

30 which are of a single frequency band, are weighted with the weight factors, such that the values of the 
weight factors may be varied for the transform factor signals of the different frequency bands. The weighted 
transform factor signals, which are of a single frequency band, are then added to one another, and an 
addition transform factor signal is thereby obtained for each of the different frequency bands. An inverse 
sub-band transfomn is then carried out on the addition transform factor signal. In this manner, as in the 

35 wavelet transform, an addition signal Sadd can be obtained from the weighted addition. 

In the embodiment described above, with the wavelet transform or the sub-band transform, the image 
signal is transformed into a plurality of image signals, each of which Is of one of a plurality of different 
frequency bands. As another alternative, the image signal may be transformed Into a plurality of image 
signals, each of which is of one of a plurality of different frequency bands, by carrying out a Fourier 

40 transform. However, for example, as illustrated in Figure 12, long filters for different frequency bands must 
be employed in the Fourier transform. Therefore, even if a fast Fourier transform process is employed, the 
apparatus for canning out the superposition processing method for a radiation image in accordance with the 
present invention will become comparatively complicated. On the other hand, in cases where the wavelet 
transform or the sub-band transform, wherein the image signal can be decomposed with a short filter, is 

45 employed, the apparatus for carrying out the superposition processing method for a radiation image in 
accordance with the present invention can be kept simple. 

The frequency dependency characteristics of the image will also vary for different doses of radiation 
delivered to the stimulable phosphor sheet. Therefore, the middle dose of radiation may be detected from 
the image read-out apparatus during the image read-out operation. Reference may be made to addition 

50 ratio tables for different doses of radiation, and the addition ratio between the Image signals may be 
determined for each of different frequency bands. 

The transfonm factor signals, which have been obtained from the wavelet transform or the sub-band 
transform described above, constitute the image signals contracted from the original image signal. For 
example, as illustrated in Figure 13. In cases where the image signal represents an image of the chest of a 

55 human body, the transform factor signals obtained from the wavelet transform or the sub-band transfomn of 
the image signal represent the images contracted from the original image. In such a chest image, the 
portions, which are to be used and therefore are required to have an appropriate image density in the 
reproduced Image, are the lung field patterns. (In Figure 13, lung field patterns 30, 30 are shown in the 
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image represented by the transform factor signal WWuv) To the portions of the stimulable phosphor sheet, 
at which the lung field patterns 30. 30 are formed, a large dose of radiation reached during the operation for 
recording the radiation image. 

As described above, the frequency dependency characteristics of the image vary for different doses of 

5 radiation delivered to the stimulable phosphor sheet. Therefore, the dose of radiation, which reached to 
each of portions of the radiation image during the operation for recording the radiation image, can be 
detected by carrying out, for example, an analysis of the probability density function of each of the 
transform factor signals, which have been obtained from the wavelet transform or the sub-band transform. 
Thereafter, the value of the weight factor with respect to the signal components of each transform factor 

10 signal representing a portion of the radiation image, to which a large dose of radiation reached during the 
operation for recording the radiation image, can be rendered larger than the value of the weight factor with 
respect to the signal components of each transform factor signal representing a portion of the radiation 
Image, to which a small dose of radiation reached during the operation for recording the radiation image. In 
this manner, an image can be obtained which has good image quality and can serve as an effective toot in, 

15 particularly, the efficient and accurate diagnosis of an illness. 

Specifically, the dose of radiation is calculated for each of the signal components of each transform 
factor signal. As for the signal components of each transform factor signal representing the picture elements 
in the lung field patterns, to which a large dose of radiation reached, the value of the weight factor is set to 
be large. For example, as for the image represented by the transfonm factor signal WWyi in Figure 13. the 

20 superposition ratio of the signal components of the transform factor signal WWui, which represent the lung 
field patterns 30, 30, to the signal components of the transform factor signal WWu, which represent the 
lung field patterns 30, 30, is set to be 2:1 . Also, the superposition ratio of the signal components of the 
transform factor signal WWui, which represent a portion 31 other than the lung field patterns 30. 30, to the 
signal components of the transform factor signal WWli, which represent the portion 31, is set to be 1:1. In 

25 this manner, the value of the weight factor with respect to the lung field patterns 30, 30, which are to be 
used and therefore are required to have an appropriate image density in the reproduced image, and to 
which a large dose of radiation reached during the operation for recording the radiation image, may be 
weighted even further, and the addition of the transform factor signals may thereby be carried out. 

In the manner described above, the value of the weight factor during the addition may also be changed 

30 in accordance with the portions of the image. In such cases, an addition signal can be obtained such that 
the value of the weight factor may be varied for different doses of radiation, which reached to portions of 
the radiation image during the operation for recording the radiation image. Therefore, an image can be 
obtained which has better image quality and can serve as a more effective tool in. particulariy, the efficient 
and accurate diagnosis of an illness. 

35 In the embodiment described above, as illustrated in Rgure 1 , the radiation images are recorded on the 
two stimulable phosphor sheets 4A and 4B. The image signals obtained from the two stimulable phosphor 
sheets 4A and 4B are then added to each other. Alternatively, as illustrated in Figure 14, the radiation image 
of the object 1 may be recorded on the single stimulable phosphor sheet 4A. As illustrated in Figure 15, two 
image signals to be added to each other may be obtained by detecting light emitted by the two surfaces of 

40 the stimulable phosphor sheet 4A, How the operation for detecting the light emitted by the two surfaces of 
the stimulable phosphor sheet 4A is carried out will be described hereinbelow. 

With reference to Figure 15, the stimulable phosphor sheet 4A is placed on endless belts 9a and 9b. 
The endless belts 9a and 9b are rotated by motors (not shown). A laser beam source 10 and a scanning 
mirror 12 are located above the stimulable phosphor sheet 4A. The laser beam source 10 produces a laser 

45 beam 11 serving as stimulating rays, which cause the stimulable phosphor sheet 4A to emit light in 
proportion to the amount of energy stored thereon during its exposure to the radiation. The scanning mirror 
12 reflects and deflects the laser beam 11, which has been produced by the laser beam source 10, such 
that the laser beam 1 1 may scan the stimulable phosphor sheet 4A in main scanning directions. A light 
guide member 14a is located above and close to the position on the stimulable phosphor sheet 4A, which is 

50 being scanned with the laser beam 11. The light guide member 14a collects light, which is emitted by the 
stimulable phosphor sheet 4A when it is scanned with the laser beam 11, from above the stimulable 
phosphor sheet 4A. Also, a light guide member 14b is located below the position on the stimulable 
phosphor sheet 4A, which is being scanned with the laser beam 11. The light guide member 14b is located 
perpendiculariy to the stimulable phosphor sheet 4A and collects the light, which is emitted by the 

55 stimulable phosphor sheet 4A when it is scanned with the laser beam 11, from below the stimulable 
phosphor sheet 4A. The light guide members 14a and 14b are located such that they may respectively be 
in close contact with photomuttipliers 15a and 15b, which photoelectrically detects the light emitted by the 
stimulable phosphor sheet 4A. The photomultipliers 15a and 15b are respectively connected to logarithmic 
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amplifiers 16a and 16b. The logarithmic amplifiers 16a and 16b are respectively connected to analog-to- 
digital converters 17a and 17b. The analog-to-digital converters 17a and 17b are connected to a memory 
18. 

The stimulable phosphor sheet 4A. on which the radiation image has been stored, is set at a 
predetermined position on the endless belts 9a and 9b. The stimulable phosphor sheet 4A, which has been 
set at the predetermined position, is conveyed by the endless belts 9a and 9b in a sub-scanning direction 
indicated by the arrow Y. Also, the laser beam 11 is produced by the laser beam source 10. The laser 
beam 11, which has been produced by the laser beam source 10, is reflected and deflected by the 
scanning mirror 12. The laser beam 11, which has thus been reflected and deflected by the scanning mirror 
12, impinges upon the stimulable phosphor sheet 4A and scans it in the main scanning directions indicated 
by the double headed arrow X. The main scanning directions are approximately normal to the sub-scanning 
direction indicated by the arrow Y. When the stimulable phosphor sheet 4A is exposed to the laser beam 
11, the exposed portion of the stimulable phosphor sheet 4A emits light in proportion to the amount of 
energy stored thereon during its exposure to the radiation. The light, which is emitted upwardly by the 
stimulable phosphor sheet 4A, is represented by reference numeral 13a. The light, which is emitted 
downwardly by the stimulable phosphor sheet 4A, is represented by reference numeral 13b. The emitted 
light 13a is guided by the light guide member 14a and photoelectrically detected by the photomultiplier 15a. 
The emitted light 13a, which has entered from the input end face of the light guide member 14a into the 
light guide member 14a, is guided through repeated total reflection inside of the light guide member 14a, 
emanates from the output end face of the light guide member 14a, and is received by the photomultiplier 
15a. The amount of the emitted light 13a representing the radiation image is converted by the photomul- 
tiplier 15a into an electric signal. In the same manner as that described above, the emitted light 13b Is 
guided by the light guide member 14b and is photoelectrically detected by the photomultiplier 15b. 

The photomultiplier 15a generates an analog output signal SA. The analog output signal SA is 
logarithmically amplified by the logarithmic amplifier 16a and converted into a digital image signal SI by the 
analog-to-digltal converter 17a. The digital image signal Si is then fed into the memory 18. Also, the 
photomultiplier 15b generates an analog output signal SB. The analog output signal SB is logarithmically 
amplified by the logarithmic amplifier 16b and converted into a digital image signal S2 by the analog-to- 
digital converter 17b. The digital Image signal S2 is then fed into the memory 18. Each of the two image 
signals SI and S2 is subjected to the wavelet transform or the sub-band transfomn and is thereby 
decomposed into a plurality of the transform factor signals, each of which is of one of a plurality of different 
frequency bands. Thereafter, in the same manner as that in the aforesaid embodiment, the transform factor 
signals, which are of a single frequency band, are weighted with the weight factors and added to each 
other. As for the transform factor signals for the high frequency band, the value of the weight factor with 
respect to the image signal Si, which has been obtained from the front surface side of the sfimulable 
phosphor sheet 4A (i.e. the side of the stimulable phosphor sheet 4A, which was closer to the radiation 
source during the operafion for recording the radiation image), is set to be larger than the value of the 
weight factor with respect to the image signal S2, which has been obtained from the back surface side of 
the sfimulable phosphor sheet 4A (i.e. the side of the sfimulable phosphor sheet 4A, which was remote fi-om 
the radiafion source during the operation for recording the radiation image). The addition transform factor 
signal, which has thus been obtained, is then subjected to the inverse wavelet transform, the inverse sub- 
band transform, or the inverse Fourier transform. As in the aforesaid embodiment, an image having good 
image quality and containing little noise component can be reproduced from the addition signal, which has 
thus been obtained. 

In the operation for detecting the light emitted by the two surfaces of the stimulable phosphor sheet, the 
stimulable phosphor sheet 4A Is scanned with the laser beam 1 1 , which has been produced by tiie single 
laser beam source 10. Alternatively, as illustrated in Figure 16. a laser beam source 10a and a scanning 
mirror 12a may be located on the front surface side of the stimulable phosphor sheet 4A. Also, a laser 
beam source 10b and a scanning mirror 12b may be located on the back surface side of the stimulable 
phosphor sheet 4A. The two surfaces of the stimulable phosphor sheet 4A may be scanned respectively 
witii laser beams 11a and 11b, which have been produced by the laser beam sources 10a and 10b. The 
light emitted by the two surfaces of the stimulable phosphor sheet 4A may thus be detected, and two image 
signals may thereby be obtained. 

In the same manner as that in the aforesaid emt>odiment. in cases where the image signals SI and S2, 
which have been obtained from the operation for detecting the light emitted by the two surfaces of the 
stimulable phosphor sheet, are subjected to the wavelet transform or the sub-band transform, the value of 
the weight factor during the addition may be changed in accordance with the portions of the image. In such 
cases, an addition signal can be obtained such that the value of the weight factor may be varied for 
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different doses of radiation, which reached to different portions of the radiation image during the operation 
for recording the radiation image. Therefore, an image can be obtained which has better image quality and 
can serve as a more effective tool in, particularly, the efficient and accurate diagnosis of an illness. 

In the embodiment described above, the processing for rendering the value of the weight factor with 

5 respect to the frequency components, which have a low signal-to-nolse ratio, smaller than the value of the 
weight factor with respect to the frequency components, which have a high signal-to-noise ratio, in 
accordance with the frequency characteristics of each of the image signals S1 and S2 is carried out with 
the wavelet transform or the Fourier transform on the Image signals S1 and S2 having been digitized by the 
analog-to-digital converter 17 or the analog-to-digital converters 17a and 17b. Alternatively, the processing 

70 described above may be carried out on the analog output signals SA and SB, which have been obtained 
from the stimulable phosphor sheets 4A and 4B in the radiation image read-out apparatus shown in Figure 
2, or which have been obtained from the stimulable phosphor sheet 4A in the radiation image read-out 
apparatus shown in Figure 15 or Figure 16. 

How the processing described above Is carried out on the analog output signals SA and SB will be 

75 described hereinbelow. 

Figure 17 is a block diagram showing how the processing is carried out on the analog output signals SA 
and SB. As described above, the analog output signals SA and SB are obtained from the stimulable 
phosphor sheets 4A and 4B in the radiation image read-out apparatus shown in Figure 2, or are obtained 
from the stimulable phosphor sheet 4A in the radiation image read-out apparatus shown in Figure 15 or 

20 Rgure 16. As illustrated in Rgure 17, the analog output signals SA and SB are respectively fed into 
logarithmic amplifiers 16a and 16b. The analog output signals SA and SB, which have been logarithmically 
amplified by the logarithmic amplifiers 16a and 16b, are respectively fed into frequency processing circuits 
40a and 40b. In the frequency processing circuits 40a and 40b, the processing described below is carried 
out. 

25 As described above, in cases where the image signals are obtained with the aforesaid operation for 
detecting light emitted by two surfaces of a stimulable phosphor sheet, the high frequency components of 
the image signal, which has been obtained from the back surface side of the stimulable phosphor sheet (i.e. 
the side of the stimulable phosphor sheet, which was remote from the radiation source during the operation 
for recording the radiation image), contain a large amount of noise components due to scattered radiation, 

30 or the like. Also, in cases where the image signals are obtained by recording the radiation images 
respectively on a plurality of stimulable phosphor sheets placed one upon another, the high frequency 
components of the image signal having been obtained from a stimulable phosphor sheet, which was remote 
from the radiation source during the operation for recording the radiation images, contain a large amount of 
noise components due to scattered radiation, or the like. Therefore, the value of the weight factor with 

35 respect to the high frequency components of the image signal, which has been obtained from the side of 
the stimulable phosphor sheet remote from the radiation source or which has been obtained from the 
stimulable phosphor sheet remote from the radiation source, is rendered smaller than the value of the 
weight factor with respect to the high frequency components of the image signal, which has been obtained 
from the side of the stimulable phosphor sheet close to the radiation source or which has been obtained 

40 from the stimulable phosphor sheet close to the radiation source. In this manner, a superposition image 
containing little noise component can be obtained. Therefore, in the frequency processing circuit 40a, the 
filtering processing for emphasizing the high frequency components of the output signal SA is carried out 
with a filter, which is shown in Rgure 18. Also, in the frequency processing circuit 40b, the filtering 
processing for reducing the high frequency components of the output signal SB is carried out with a filter, 

45 which Is shown in Figure 19. Because the output signals SA and SB are the analog signals, the filtering 
processing described above is carried out only on the signal components of each of the output signals SA 
and SB with respect to the main scanning direction. 

The output signals SA and SB, which have been subjected to the frequency processing in the 
frequency processing circuits 40a and 40b, are respectively fed into aliasing eliminating filters 41a and 41b. 

50 In the aliasing eliminating filters 41 a and 41 b, errors due to aliasing are removed from the output signals SA 
and SB. Thereafter, the output signals SA and SB are fed into analog-to-digital converters 17a and 17b and 
converted into digital image signals SI and S2. The digital image signals SI and S2, which have thus been 
obtained, are added to each other. Also, in the same manner as that in the aforesaid embodiment, 
predetermined image processing is carried on the addition signal Sadd In the image processing means. The 

55 addition signal Sadd, which has been obtained from the image processing, is then used for reproducing a 
visible image in the image reproducing means. 

The filters shown in Rgures 18 and 19 are set such that the DQE of the addition signal Sadd resulting 
from the addition of the digital Image signals Si and S2 may become largest. Specifically, the filters shown 
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in Figures 18 and 19 are set such that the addition ratio at 1 cycle/mnr) may be SA:SB = 0.5:0.5, the addition 
ratio at 2 cycles/mnr^ may l^e SA:SB = 0.7:0.3, and the addition ratio at 3 cycles/mm may be SA:SB = 0.9:0.1. 
This means that the response of the means, such as the frequency processing circuit 40b and the aliasing 
eliminating filter 41b, for carrying out the filtering processing on the output signal SB. with respect to the 
high frequency components of the output signal SB at approximately 3 cycles/mm may be approximately 
20% of the response with respect to the low frequency components of the output signal SB at 1 cycle/mm. 
Therefore, it is sufficient for the means, such as the frequency processing circuit 40b and the aliasing 
eliminating filter 41b, for carrying out the filtering processing on the output signal SB. to be constituted of 
circuits having a narrower processing range than in the means, which are ordinarily employed in 
apparatuses for carrying out the image read-out operations. 

In the manner described above, the processing for rendering the value of the weight factor with respect 
to the frequency components, which have a low signal-to-noise ratio, smaller thcin the value of the weight 
factor with respect to the frequency components, which have a high signal-to-noise ratio, In accordance with 
the frequency characteristics of each of the output signals SA and SB is carried out on the analog output 
signals SA and SB. As a result, as in the cases wherein the processing described above is carried out on 
the digital image signals SI and S2, the addition signal having a high signal-to-noise ratio over the entire 
frequency bands can be obtained. Therefore, a superposition image having good image quality can be 
reproduced from the addition signal. Also, the high frequency components of the output signal SB are 
reduced sufficiently, and therefore the cut-off frequency of the aliasing eliminating filter 41 b can be set to 
be low. As a result, the aliasing can be reduced sufficiently during the analog-to-digital conversion carried 
out in the analog-to-digital converter 17b. Further, the parts constituting the logarithmic amplifier 16b and 
the aliasing eliminating filter 41b for the output signal SB are not required to carry out quick processing. 
Therefore, these parts may be constituted of cheap parts, such as operational amplifiers and transistors, 
which cannot achieve quick processing. Accordingly, the cost of the apparatus can be kept low. 

In the aforesaid operation for carrying out the processing on the analog output signals SA and SB, the 
frequency processing is carried out by the frequency processing circuits 40a and 40b. Alternatively, the 
effects of the frequency processing circuits 40a and 40b for emphasizing the high frequency components of 
the output signal SA and reducing the high frequency components of the output signal SB may be achieved 
by changing the frequency characteristics of the logarithmic amplifiers 16a and 16b and the aliasing 
eliminating filters 41a and 41b. 

Furthermore, in the aforesaid operation for carrying out the processing on the analog output signals SA 
and SB, the processing Is carried out on both of the analog output signals SA and SB. Alternatively, 
frequency processing may be carried out on either one of the analog output signals SA and SB. The 
processed analog output signal and the unprocessed analog output signal may be converted into digital 
image signals, and the digital image signals may then be added to each other. 

In the embodiment described above, when the two image signals are to be superposed one upon the 
other, each of the image signals are decomposed into a plurality of the transform factor signals, each of 
which is of one of a plurality of different frequency bands. The transform factor signals are then weighted. 
Alternatively, the weighting may be carried out in the manner described above when energy subtraction 
processing is carried out on two image signals. How the image signals to be subjected to the energy 
subtraction processing are weighted will be described hereinbelow. 

Rgure 20 shows an image recording apparatus for carrying out the so-called "one-shot energy 
subtraction processing", in which radiation 2 having passed through a single object 1 is delivered to a first 
stimulable phosphor sheet 4A, the energy distribution of the radiation 2 after having passed through the first 
stimulable phosphor sheet 4A is changed, and the radiation 2 now having the different energy distribution Is 
delivered to a second stimulable phosphor sheet 4B. Specifically, the first stimulable phosphor sheet 4A is 
located at the position close to a radiation source 3, and the second stimulable phosphor sheet 4B is 
located at the position slightly spaced from the first stimulable phosphor sheet 4A. A radiation energy 
converting filter 5, which may be constituted of a copper plate, is located between the two stimulable 
phosphor sheets 4A and 4B. In this state, the radiation source 3 is driven. In this manner, on the first 
stimulable phosphor sheet 4A. a radiation image of the object 1 is recorded with the radiation 2 containing 
soft radiation. Also, on the second stimulable phosphor sheet 4B, a radiation image of the object 1 is 
recorded with the radiation 2 free from soft radiation. 

In the manner described above, the two radiation images, In which different Images of a specific 
structure of the object 1 are embedded, are stored on the two stimulable phosphor sheets 4A and 4B. 

Thereafter, the radiation images are respectively read out from the two stimulable phosphor sheets 4A 
and 4B by using the radiation image read-out apparatus shown in Figure 2. Two digital image signals Si 
and S2 representing the radiation images are thereby obtained. The obtained Image signals SI and S2 are 
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stored in the storage medium 18. 

Subtraction processing Is then carried out on the digital image signals SI and S2, which have been 
obtained in the manner described above. Figure 21 is a block diagram showing an apparatus for carrying 
out a first embodiment of the energy subtraction processing method in accordance with the present 

5 invention. First, the image signals SI and 82 are read from an image file ISA and an image file 18B in the 
storage medium 1 8 and fed into a wavelet transform means 30. The wavelet transfonm means 30 carries out 
a wavelet transform of each of the two image signals SI and S2. In the same manner as that in the 
aforesaid embodiment of the superposition processing method for a radiation image in accordance with the 
present invention, each of the two image signals SI and S2 is thus decomposed into a plurality of wavelet 

70 transform factor signals, each of which is of one of a plurality of different frequency bands. In this manner, 
wavelet transform factor signals WWuo through WWun, WVuo through WVun. VWuo through VWun. and 
WuN are obtained from the image signal S1 . Also, wavelet transfonm factor signals WWlo through WWln. 
WVlo through WVln, \A/Vlo through VW^n, and Wln are obtained from the image signal S2. 

The wavelet transform factor signals, which have thus been obtained, are fed into a weighting and 

75 subtraction means 31. In the weighting and subtraction means 31, subtraction processing Is carried out 
such that the value of the weight factor with respect to the frequency band, which contains much noise 
components, may be rendered smaller than the value of the weight factor with respect to the frequency 
band, which contains less noise components. The information in the high frequency band of the image 
signal S2 contains noise, due to scattered radiation during the image recording operation, or the like. Also, 

20 the information on the high frequency band side of the image signal S2 has been rendered unsharp due to 
the location of the stimulable phosphor sheet 4B remote from the radiation source. Therefore, the two image 
signals Si and S2 have the frequency characteristics shown in Figures 6A and 6B. Accordingly, when the 
subtraction processing is carried out on the two image signals SI and S2, the values of the weight factors 
are determined in accordance with the weight table shown in Figure 9. The wavelet transform factor signals, 

25 which are of a single frequency band, are weighted in accordance with the weight table shown in Figure 9, 
and the weighted wavelet transform factor signals, which are of a single frequency band, are then 
subtracted from one another. 

Specifically, the weighted subtractions of the wavelet transform factor signals, which are of a single 
frequency band, are carried out with Formula (13) 

30 

WW, = t • WWut - (1 ■ t) WWli (13) 
WV| = t • WVui - (1 - 1) WVu 
VW, = t • VWui - (1 • t) \AA/u 
W, = t • Wui - (1 - 1) Wui 

35 

In the manner described above, subtraction wavelet transform factor signals WWi through WWn, WVi 
through WVn, VWi through VW^, and Wi through Wn are obtained in the weighting and subtraction 
means 31. Thereafter, as in the aforesaid embodiment of the superposition processing method for a 
radiation image In accordance with the present invention, in an inverse wavelet transform means 32, an 

40 inverse wavelet transform is carried out on each of the subtraction wavelet transform factor signals. A 
subtraction signal, which has been obtained from the inverse wavelet transform, is subjected to predeter- 
mined image processing in an image processing means 33, fed into an image reproducing means 34, and 
used for the reproduction of a visible radiation image. 

In the manner described above, the wavelet transform is carried out on each of the two image signals 

45 SI and S2, and each of the two image signals SI and S2 is thereby decomposed into a plurality of the 
wavelet transform factor signals, each of which is of one of a plurality of different frequency bands. The 
wavelet transform factor signals, which are of a single frequency band, are weighted with the weight factors, 
such that the values of the weight factors may be varied for the wavelet transfonn factor signals of the 
different frequency bands. As for the wavelet transform factor signals for the high frequency band, the value 

50 of the weight factor with respect to the image signal, which has been obtained from the stimulable phosphor 
sheet located at the position closer to the radiation source, is set to be larger than the value of the weight 
factor with respect to the image signal, which has been obtained from the stimulable phosphor sheet 
located at the position remote from the radiation source. In this manner, a subtraction signal containing little 
noise can be obtained. From the subtraction signal thus obtained, an image having good image quality and 

55 containing little noise component can be reproduced. 

In the aforesaid embodiment of the energy subtraction processing method in accordance with the 
present invention, the so-called "one-shot energy subtraction processing" is carried out. In which the two 
radiation images are recorded on the stimulable phosphor sheets with a single, simultaneous exposure to 
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the radiation, and the two image signals SI and S2 are obtained from the stimulable phosphor sheets. 
Alternatively, two image signals may be obtained with the so-called "two-shot energy subtraction process- 
ing." in which two stimulable phosphor sheets are exposed one after another to two kinds of radiation 
having different energy distributions. In cases where the subtraction processing is carried out on the two 

6 image signals obtained with the two-shot energy subtraction processing, each of the two image signals 
obtained may be decomposed into a plurality of the transform factor signals, each of which is of one of a 
plurality of different frequency bands. The transform factor signals, which are of a single frequency band, 
may be weighted and subtracted from each other. In such cases, the values of the weight factors may be 
determined in accordance with, for example, the MTPs which are obtained from the energy distributions of 

10 the two kinds of the radiation employed when the radiation images were recorded on the two stimulable 
phosphor sheets. 

In the aforesaid embodiment of the energy subtraction processing method in accordance with the 
present invention, each of the two Image signals is decomposed by the wavelet transform into a plurality of 
the wavelet transform factor signals, each of which is of one of a plurality of different frequency bands. 

75 Alternatively, each of the two image signals may be decomposed by the sub-band transfomn or the Fourier 
transform into a plurality of the transform factor signals, each of which is of one of a plurality of different 
frequency bands. However, for example, as illustrated in Figure 12, long filters for different frequency bands 
must be employed in the Fourier transform. On the other hand, in cases where the wavelet transform or the 
sub-band transform, wherein the image signal can be decomposed with a short filter, is employed, the 

20 apparatus for carrying out the energy subtraction processing method in accordance with the present 
invention can be kept simple. 

The transform factor signals, which have been obtained from the wavelet transform or the sub-band 
transform described above, constitute the image signals contracted from the original image signal. For 
example, as illustrated in Figure 13, in cases where the image signal represents an image of the chest of a 

25 human body, the transform factor signals obtained from the wavelet transform or the sub-band transform of 
the image signal represent the images contracted from the original image. Therefore, as in the aforesaid 
embodiment of the superposition processing method for a radiation image in accordance with the present 
invention, the value of the weight factor during the subtraction processing may also be changed in 
accordance with the portions of the image. In such cases, a subtraction signal can be obtained such that the 

30 value of the weight factor may be varied for different doses of radiation, which reached to portions of the 
radiation image during the operation for recording the radiation image. Therefore, an image can be obtained 
which has better image quality and can serve as a more effective tool in. particularly, the efficient and 
accurate diagnosis of an illness. 

In the aforesaid embodiment of the energy subtraction processing method in accordance with the 

35 present invention, as illustrated in Rgure 17, the analog output signals SA and SB may be subjected to the 
frequency processing with the filters shown in Figures 18 and 19 and may then be converted into digital 
image signals SI and S2. The subtraction processing may then be carried out on the digital image signals 
SI and S2. 

Also, in the aforesaid operation for carrying out the processing on the analog output signals SA and SB, 

40 the processing is carried out on both of the analog output signals SA and SB. Alternatively, frequency 
processing may be carried out on either one of the analog output signals SA and SB. The processed analog 
output signal and the unprocessed analog output signal may be converted into digital image signals, and 
the digital image signals may then be subtracted from each other. 

A second embodiment of the superposition processing method for a radiation image in accordance with 

45 the present invention and a second embodiment of the energy subtraction processing method in accor- 
dance with the present invention will be described hereinbelow. 

Image signals S1 and S2 are obtained in the same manner as that described above with reference to 
Figures 1 and 2. Thereafter, a superposition processing is carried out on the image signals SI and S2. 
Rgure 25 is a block diagram showing an apparatus for carrying out the second embodiment of the 

50 superposition processing method for a radiation image in accordance with the present invention. First, the 
image signals SI and S2 are read from an image file 18A and an image file IBB in the storage medium 18 
and fed into an image processing means 119. The image processing described below is carried out on the 
two image signals SI and S2, which have been fed into the image processing means 119. The image 
signals, which have been obtained from the image processing, are fed into a weighting and superposition 

55 means 120 and subjected to a weighted addition. An addition signal, which has been obtained from the 
weighting and superposition means 120, is fed into an image reproducing means 121, such as a CRT 
display device, and is used for reproducing a visible image. 
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How the image processing is carried out in the image processing means 119 will be described 
hereinbetow. 

As described above with reference to Rgures 6A through 9. the optimum addition ratio, which gives the 
largest DQE serving as an index of an image, varies for different frequencies. Also, as described above, the 
image signal S1 carries the information up to the high frequency band. However, in the image signal S2, 
noise components are predominant in the information of the high frequency band. Therefore, as for the 
image signal SI , the response with respect to the high frequency band is emphasized. As for the image 
signal S2. the response with respect to the high frequency band is reduced. The image signals SI and S2 
are then added to each other. In this manner, an addition signal can be obtained which represents a 
radiation image having good image quality. 

How the image signals 81 and S2 are processed will be described hereinbelow. 

First, an unsharp mask signal Susi for the image signal SI is calculated. How the unsharp mask signal 
Susi is calculated will be described hereinbelow. 

Figure 26 is an explanatory view showing picture elements in an image and image signal components 
of the image signal SI, which represent the picture elements. In Figure 26, dots represent the picture 
elements, and the symbols, such as SiJ, represent the image signal components of the image signal Si, 
which represent the con'esponding picture elements. 

The value of an unsharp mask signal Si| for the picture element, which is located In the middle of the 
region surrounded by the chained line, is calculated with the fonmula 



I J 

Susi =■ 



2 



£ S I + k 



J + I 



M • N 



(14) 



By can7ing out the calculation for every picture element, the unsharp mask signal Susi for the whole image 
is generated. The values of M and N are determined arbitrarily in accordance with the intervals, with which 
the image signal SI is sampled, the characteristics of the radiation image, the type of the desired image 
processing, or the like. 

After the unsharp mask signal Susi is obtained, the process shown below is carried out on the unsharp 
mask signal Susi . 

SV = SI + ^ (SI - Susi) (15) 

where 0 represents the degree of emphasis. 

Specifically, the process represented by Formula (15) is carried out on the image signal SI by taking 
the value of as the value for emphasizing the high frequency band of the image signal SI . In this manner, 
the high frequency band of the image signal SI is emphasized. 

On the other hand, as for the image signal S2, an unsharp mask signal Sus2 for the image signal S2 is 
calculated with Formula (14). The unsharp mask signal Sus2 Is taken as a signal S2\ which Is to be added 
to the image signal SI. Specifically, as for the image signal S2, the image is made unsharp by carrying out 
the unsharp mask processing with Formula (14). In this manner, the response with respect to the high 
frequency band of the image signal S2 is reduced. 

The image signals SI and S2 are processed in the manner described above, and image signals SV 
and S2' are thereby obtained. The image signals SI' and S2' are weighted with predetermined values of the 
weight factors and are thereafter added to each other. Specifically, an addition signal Sadd is obtained with 
the formula 

Sadd = t • Sr + (1 - 1) • S2* (16) 

The addition signal Sadd obtained in this manner is fed into an image reproducing means 121 and used 
for reproducing a visible Image. 

The image reproducing means may be a display device, such as a cathode ray tube (CRT) display 
device, or a recording apparatus for recording an image on photosensitive film by scanning the photosen- 
sitive film with a light beam. 

In the manner described above, the frequency emphasis processing is carried out on the image signal 
Si, and the Image signal SV is thereby obtained. Also, the unsharp mask processing is carried out on the 
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image signal S2, and the image signal S2* is thereby obtained. Thereafter, the innage signals SV and S2* 
are added to each other. This operation is substantially equivalent to the operation, in which the addition is 
carried out such that the addition ratio may be different between the low frequency band and the high 
frequency band of the image signals. In this manner, the addition signal can be obtained such that the noise 

5 components obtained from the stimulable phosphor sheet, which was located at the position remote from 
the radiation source, may be reduced, and such that the information of the high frequency band obtained 
from the stimulable phosphor sheet, which was located at the position close to the radiation source, may be 
emphasized. Therefore, from the addition signal thus obtained, an image having good image quality and 
containing tittle noise component can be reproduced. Also, the amount of calculation is not so large as in 

70 the wavelet transform or the Fourier transform. Therefore, the apparatus for carrying out the second 
embodiment of the superposition processing method for a radiation image in accordance with the present 
invention can be kept simple. Further, the operation can be carried out quickly. 

In the aforesaid second embodiment of the superposition processing method for a radiation image in 
accordance with the present invention, the image processing is earned out on both of the Image signals Si 

75 and S2. Alternatively, the image processing described above may be carried out on either one of the image 
signals S1 and S2. However, in cases where the image processing is carried out on both of the image 
signals S1 and S2, an addition signal representing an image having better image quality can be obtained. 

Also, in the aforesaid second embodiment of the superposition processing method for a radiation image 
in accordance with the present invention, the unsharp mask signal is calculated with Formula (14) by using 

20 the mean value of the values of the image signal components representing the picture elements located 
within the mask having a predetermined range. Alternatively, a weighted mean value of the values of the 
image signal components representing the picture elements located within the mask may be employed 
such that the frequency characteristics of the image signal can be controlled. As another alternative, any of 
values, which are representative of the values of the image signal components representing the picture 

25 elements located within the mask, such as the median in the mask, may be employed. 

Further, in the aforesaid second embodiment of the superposition processing method for a radiation 
image in accordance with the present invention, as the image processing of the image signals Si and S2, 
the frequency emphasis processing utilizing the unsharp mask signal and the unsharp mask processing are 
carried out. Alternatively, other types of frequency processing may be carried out. How the frequency 

30 processing is carried out on the image signals SI and S2 in a modification of the second embodiment of 
the superposition processing method for a radiation image in accordance with the present invention will be 
described herelnbelow. 

If the image processing carried out on the image signal SI is expressed as F1(S1), and the image 
processing carried out on the image signal S2 is expressed as F2{S2). the addition signal Sadd in the 
35 aforesaid second embodiment of the superposition processing method for a radiation image in accordance 
with the present invention can be expressed as 

Sadd = t • F1(S1) + (1 - 1) • F2(S2) (17) 

40 Specifically, F1{S1) represents the frequency emphasis processing of the image signal SI, and F2{S2) 
represents the unsharp mask processing of the image signal S2. In lieu of the processing F1(S1) and the 
processing F2{S2), the convolution of each of the image signals SI and S2 may be carried out with masks 
having predetermined frequency characteristics. In this manner, signal-to-noise ratio of the addition signal 
can be kept high. 

45 Specifically, the convolution of the image signal S2 is carried out with a filter F2 having the 
characteristics shown in Figure 27. Also, the convolution of the image signal SI is carried out with a filter F1 
having the characteristics shown in Figure 28. In this manner, the response with respect to the high 
frequency band of the image signal S2 can be restricted, and the response with respect to the high 
frequency band of the image signal SI can be emphasized. Therefore, an addition signal representing an 

50 image having good image quality can be obtained. 

The filter F2 having the characteristics shown in Figure 27 is a one-dimensional tilter and has the 
coefficients 

-2/122 -1/122 17/122 31/122 17/122 -1/122 -2/122 

If the one-dimensional filter is expressed as A(i), and a two-dimensional filter is expressed as Aij, the 
55 following formula can obtain with respect to the image. 

Aij = A(i) X A(i) / EA(n) (18) 
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Specifically, the convolution of the image signal is carried out with the one-dimensional filter A(i) with 
respect to the X direction, and thereafter the convolution is carried out with respect to the Y direction. 
Thereafter, the value obtained from the convolution is divided by the total sum of the one-dimensional filter. 
In this manner, the two-dimensional filter can be obtained. The values of the two-dimensional filter are 
5 shown below (at this time, the total sum is 59/1 22). 

Table 4 



10 



15 



4/7198 


2/7198 


-34/7198 


-62/7198 


-34/7198 


2/7198 


4/7198 


2/7198 


1/7198 


-17/7198 


-31/7198 


-17/7198 


1/7198 


2/7198 


-34/7198 


-17/7198 


289/7198 


527/7198 


289/7198 


-17/7198 


-34/7198 


-62/7198 


-31/7198 


527/7198 


961/7198 


,527/7198 


-31/7198 


-62/7198 


-34/7198 


-17/7198 


289/7198 


527/7198 


289/7198 


-17/7198 


-34/7198 


2/7198 


1/7198 


-17/7198 


-31/7198 


-17/7198 


1/7198 


2/7198 


4/7198 


2/7198 


-34/7198 


-62/7198 


-34/7198 


2/7198 


4/7198 



Also, as shown in Formula (17), in the aforesaid second embodiment of the superposition processing 
method for a radiation image in accordance with the present invention, the image signals obtained from the 

20 processing are respectively multiplied by t and 1-t such that the energy of the signal values of the addition 
signal may become Identical with the energy of the original image signals SI and S2. In this modification of 
the second embodiment of the superposition processing method for a radiation image in accordance with 
the present invention, the values of the filters Fl and F2 are set appropriately such that the image signals 
obtained from the processing are may not be multiplied by t and 1-t as shown in Formula (17). Specifically, 

25 in this modification, the filters Fl and F2 are set such that Fl + F2 = 1. The values of the filter F2 are set as 
described above. How the values of the filter Fl are set such that Fl + F2 = 1 will be described hereinbelow. 

First, the median (31/122) of the filter F2 is altered to a value of (1 - median), and the signs of the other 
filter elements are reversed. In this manner, the filter Fl having the values shown below can be obtained. 
2/122 1/122 -17/122 91/122 -17/122 1/122 2/122 The filter F1 has the characteristics shown in Figure 28. 

30 Also, as for the two-dimensional filter, in the same manner as that described above, the median of the 
filter is altered to a value of (1 - median), and the signs of the other filter elements are reversed. In this 
manner, a two-dimensional filter Fl having the values shown below can be obtained. 

Table 5 
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40 



-4y7198 


-2/7198 


34/7198 


62/7198 


34/7198 


-2/7198 


-4/7198 


-2>7198 


-1/7198 


17/7198 


31/7198 


17/7198 


-1/7198 


-2/7198 


34/7198 


17/7198 


-289/7198 


-527/7198 
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45 After the filters Fl and F2 have been obtained in the manner described above, the convolution of the 
image signal Si with the filter Fl and the convolution of the image signal S2 with the filter F2 are carried 
out. Thereafter, the addition is carried out. In this manner, as In the aforesaid second embodiment of the 
superposition processing method for a radiation image in accordance with the present invention, an image 
signal can be obtained in which the response with respect to the high frequency band of the image signal 

50 81 has been emphasized. The processing described above can be represented by Formula (19). 

Sadd = F, - Si + F2 * S2 (19) 

where FTSI represents the convolution of SI with Fl. 
55 The convolution of the image signal SI with the filter Fl and the convolution of the image signal 82 with 
the filter F2 are thus carried out, and the obtained values are then added to each other. In this manner, the 
addition with respect to each of different frequencies can be carried out such that 
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as for 0 cycle 


S1:S2 


= 0.516:0.484 


as for 1 cycle 


S1:S2 


= 0.515:0.485 


as for 2 cycles 


S1:S2 


= 0.620:0.380 


as for 3 cycles 


S1:S2 


= 0.845:0.155 


as for 4 cycles 


S1:S2 


= 0.987:0.013 



In cases where the frequency characteristics of each of the image signals SI and S2 vary for the X 
direction and the Y direction, a different filter may be employed for each direction. In such cases, if the 
10 different filter is expressed as B(j), the two-dimensional filter can be represented by Formula (20). 

Aij = A(j) X B(j) / EA(n) (20) 

However, in such cases, it is necessary that i:A(n) = EB(n). 
/5 For reasons of the apparatus, for example, in cases where the capacity of the memory in the apparatus 
is small, it may be desired to carry out the convolution of the image signal Si and the convolution of the 
image signal S2 with a single filter. In such cases, the addition may be can'ied out with the formulas shown 
below. 

20 Sadd = (Si - F2 * Si ) + (Fa * S2) (21) 

(Because F1 + F2 = 1 . , 

Sadd = Fi * Si + F2 * S2 (22) 
25 = (1 - F2) • Si + F2 * S2 

= (Si ■F2*Sl) + (F2-S2)) 
or Sadd = Fi * Si + S2 - Fi ' S2 

In this manner, it is sufficient for only a single kind of filter to be stored in the memory. In such cases, the 
30 same results as those with Formula (19) can be obtained. 

Further, in such cases. Formulas (21 ) and (22) can be rewritten as 

Sadd = Si + F2 *(S2 - Si) (2r) 
Sadd = S2 + Fi *{Si -S2) (22') 

35 

Specifically, the same results as those with Formula (19) can be obtained by calculating the difference 
between the image signals SI and S2 and carrying out the convolution of the difference value with the filter 
F1 or the filter F2. In such cases, as in the operations with Formulas (21) and (22), it is sufficient for only a 
single kind of filter to be stored in the memory. Furthermore, it is sufficient for the convolution of the 
40 difference value between the image signals S1 and S2 to be canied out only one time with the filter F1 or 
the filter F2. Therefore, the apparatus can be kept simple, Its cost can be kept low. and the operation can 
be carried out quickly. 

In the aforesaid modification of the second embodiment of the superposition processing method for a 

radiation image in accordance with the present invention, in order for the image signals S1 and S2 to be 
45 added to each other in a desired addition ratio, it is necessary for the value of the filter coefficient to be 

determined such that the frequency characteristics of the filter F2 may become the desired characteristics. 

How the filter F2 is determined will be described hereinbelow. 

As a first technique for determining the filter F2, the Fourier transform Is carried out on an arbitrary 

value of the filter coefficient, and the frequency characteristics are investigated. The value of the filter 
50 coefficient is then altered, the Fourier transform is earned out on the new value of the filter coefficient, and 

the frequency characteristics are investigated. The value of the filter coefficient is thus finely adjusted, trial 

and error are repeated, and the filter is thereby determined such that it may have the desired frequency 

characteristics. 

As a second technique for determining the filter F2, equations are created from sets of desired 
55 frequencies and the responses by taking the filter coefficient as an undetermined value, and simultaneous 
linear equations are thereby obtained. The simultaneous linear equations 
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00 

R(fo 2^_^a (n) • exp (-j • 2;r*fo -n .T)...(23) 

5 can be obtained, where fo represents the frequency. R(fo) represents the desired response, a(n) represents 
the fitter coefficient, and T represents the sampling interval. The approximate solution of the filter coefficient 
a(n) satisfying the simuttaneous linear equations is determined with the method of least squares. 

In the aforesaid second embodiment of the superposition processing method for a radiation image In 
accordance with the present invention, as illustrated in Figure 1, the radiation images are recorded on the 

70 two stimulable phosphor sheets 4A and 4B. The image signals obtained from the two stimulable phosphor 
sheets 4A and 4B are then added to each other. Alternatively, as Illustrated in Figure 14, the radiation image 
of the object 1 may be recorded on the single stimulable phosphor sheet 4A. As illustrated in Figure 15, two 
image signals to be added to each other may be obtained by detecting light emitted by the two surfaces of 
the stimulable phosphor sheet 4A. 

;5 The two image signals SI and S2, which have been obtained from the analog-to-digital converters 17a 
and 17b in the same manner as that described above with reference to Figure 15, are added to each other 
in the same manner as that in the aforesaid second embodiment of the superposition processing method for 
a radiation image in accordance with the present invention. In this manner, as in the aforesaid second 
embodiment of the superposition processing method for a radiation image in accordance with the present 

20 invention, an image having good image quality and containing little noise component can be reproduced 
from the addition signal, which has thus been obtained. 

As another alternative, as illustrated in Figure 16, the two surfaces of the stimulable phosphor sheet 4A 
may be scanned respectively with laser beams 11a and 11b, which have been produced by the laser beam 
sources 10a and 10b. The light emitted by the two surfaces of the stimulable phosphor sheet 4A may thus 

25 be detected, and two image signals may thereby be obtained. 

In cases where the operation for quickly reading out a radiation image is carried out in the aforesaid 
second embodiment of the superposition processing method for a radiation image in accordance with the 
present invention, the stimulating rays are moved very quickly on the stimulable phosphor sheet. Therefore, 
with certain kinds of stimulable phosphors constituting the stimulable phosphor sheets, the problems often 

30 occur in that the stimulable phosphor sheet cannot emit light immediately after being exposed to the 
stimulating rays, and a time lag occurs between when the stimulable phosphor sheet is exposed to the 
stimulating rays and when the stimulable phosphor sheet emits light. At an image contour portion, or the 
like, the amount of the light emitted by the stimulable phosphor sheet changes sharply. Therefore, the 
values of the image signal, which correspond to an image contour portion, or the like, should change 

35 sharply in the main scanning direction. However, if the time lag occurs between when the stimulable 
phosphor sheet is exposed to the stimulating rays and when the stimulable phosphor sheet emits light, the 
values of the obtained image signal, which correspond to an image contour portion, or the like, will not 
change sharply in the main scanning direction. As a result, the visible image reproduced from the image 
signal becomes unsharp in the main scanning direction, and the sharpness of the reproduced image cannot 

40 be kept high. 

Even in such cases, with the modification of the aforesaid second embodiment of the superposition 
processing method for a radiation image in accordance with the present invention, the sharpness of the 
reproduced image can be prevented from becoming low by changing the kinds of the filters. How the 
sharpness of the reproduced image is prevented from becoming low will be described hereinbelow. 

45 In cases where the radiation is delivered simultaneously to the two stimulable phosphor sheets 4A and 
4B placed one upon the other as shown in Figure 1 , the radiation images being thereby stored respectively 
on the stimulable phosphor sheets 4A and 48, and the two image signals SI and S2 are obtained 
respectively from the stimulable phosphor sheets 4A and 48 by using the radiation image read-out 
apparatus shown in Figure 2, the image signal SI carries the information up to the high frequency band. 

50 However, in the image signal S2, noise components are predominant in the information of the high 
frequency band. 

Also, in cases where each of the stimulable phosphor sheets 4A and 48 is quickly (e.g., at a read-out 
clock of at least 1 .2us) scanned with the stimulating rays during the image read-out operation, the change 
in the values of each of the image signals SI and S2 will become unsharp in the main scanning direction 
55 (i.e., in the X direction). As a result, the visible image reproduced from the image signal becomes unsharp 
with respect to the main scanning direction. The problems with regard to the reproduced image becoming 
unsharp with respect to the main scanning direction occur particulariy markedly with respect to the image 
signal SI . 
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Therefore, in order to prevent the reproduced inr^age fronn becoming unsharp with respect to the main 
scanning direction, the convolution of the Image signal SI with respect to the main scanning direction and 
the convolution of tiie image signal S1 with respect to tiie sub-scanning direction are carried out with 
different filters such that the high frequency components of the image signal Si carrying the information of 

5 the high frequency components can be emphasized even further with respect to the main scanning 
direction. In the image signal S2, in which noise components are predominant in the information of the high 
frequency band, the change in the values of the image signal 82 is more unsharp in both of the main 
scanning direction and the sufc>-scanning direction than in the image signal 81. Therefore, ttie convolution of 
the image signal 82 with respect to the main scanning direction and the convolution of the image signal 82 

TO with respect to tiie sub-scanning direction are carried out with filters such that the response with respect to 
the high frequency components can be restricted. 

Specifically, the convolution of the image signal 81 with respect to the main scanning direction is 
carried out witti a tilter F3 shown in Figure 29. The convolution of tfie Image signal 82 with respect to the 
main scanning direction is carried out with a filter F4 shown in Figure 30. Also, the convolution of the image 

75 signal 81 with respect to the sub-scanning direction is earned out with the filter F1 shown in Figure 28, and 
the convolution of the image signal 82 with respect to the sub-scanning direction is carried out with the tilter 
F2 shown in Figure 27. As illustrated in Figure 31, as for the high frequency band, the sum of the titters F3 
and F4 shown in Figures 29 and 30 are F3 + F4>1. Therefore, In cases where the convolution of the image 
signal 81 with respect to the main scanning direction is carried out with the tilter F3 shown in Figure 29, the 

20 convolution of the image signal 82 with respect to the main scanning direction is carried out with the tilter 
F4 shown in Figure 30, and the image signals 81 and 82 obtained from the processing are then added to 
each other, the high frequency band is emphasized markedly in the resulting addition signal with respect to 
the main scanning direction. Accordingly, even if tiie image signals 81 and 82 are unsharp in the main 
scanning direction, tiie high frequency components of the addition signal can be emphasized, and the 

25 sharpness with respect to tiie main scanning direction can be kept high. In this manner, the reproduced 
image can be prevented from becoming unsharp with respect to the main scanning direction, and the 
reproduced image having a high sharpness can be obtained. 

In cases where the main scanning is carried out quickly during the operation for detecting light emitted 
by the two surfaces of the stimulable phosphor sheet, as in the superposition processing, the sharpness of 

30 the image signal in the main scanning direction becomes low. Therefore, as for the image signals 81 and 
82 obtained from the two surfaces of the stimulable phosphor sheet, the convolution is carried out with the 
filters F3 and F4. In this manner, as in the superposition processing, the reproduced image can be 
prevented from becoming unsharp with respect to the main scanning direction, and the reproduced image 
having a high sharpness can be obtained. 

35 In the aforesaid second embodiment of the superposition processing method for a radiation image in 
accordance with the present invention, the two image signals are superposed one upon the other. The 
processing described above can also be carried out when energy subtraction processing is carried out on 
two image signals. How the image processing is carried out on the image signals to be subjected to the 
energy subtraction processing will be described hereinbelow. 

40 In the image recording apparatus for carrying out the one-shot energy subtraction processing, which is 
shown in Figure 20, in the same manner as that described above, the two radiation images, in which 
different images of a specific structure of the object 1 are embedded, are stored on the two stimulable 
phosphor sheets 4A and 4B. 

Thereafter, tiie radiation images are respectively read out from tiie two stimulable phosphor sheets 4A 

45 and 4B by using the radiation image read-out apparatus shown in Figure 2. Two digital image signals 81 
and 82 representing the radiation images are thereby obtained. The obtained image signals 81 and 82 are 
stored in the storage medium 1 8. 

Subti'action processing is then carried out on the digital image signals 81 and 82, which have been 
obtained in the manner described above. Figure 32 is a block diagram showing an apparatus for carrying 

50 out a second embodiment of the energy subtraction processing method In accordance with the present 
Invention. First, the image signals 81 and 82 are read from an image file 18A and an Image file 18B in the 
storage medium 18 and fed into an image processing means 130. In the same manner as that in the 
aforesaid second embodiment of the superposition processing method for a radiation image In accordance 
with the present Invention, the two image signals 81 and 82 are subjected to teh processing, such as the 

55 unsharp mask processing, teh frequency emphasis processing, and the convolution with the mask, in the 
image processing means 130. Processed image signals 81' and 82' are tfiereby obtained. 

The image signals SV and 82' thus obtained are fed Into a weighting and subtraction means 131. In the 
weighting and subtraction means 131, subtraction processing is carried out on the Image signals 81' and 



39 



EP 0 677 780 A2 



S2\ 

Specifically, the weighted subtractions are carried out with Formula (24) 

Ssub = tl . Sr - 12 . S2' (24) 

5 

where t1 and t2 represent the weight factors. 

In the manner described above, the subtraction signal Ssub is obtained from the weighting and 
subtraction means 131. The subtraction signal Ssub is fed into an image reproducing means 132 and used 
for reproducing a visible image. 

70 In the aforesaid second embodiment of the energy subtraction processing method in accordance with 
the present invention, the one-shot energy subtraction processing is carried out. The aforesaid second 
embodiment of the energy subtraction processing method in accordance with the present invention is also 
applicable to the two-shot energy subtraction processing. 

In the manner described above, the frequency emphasis processing, the unsharp mask processing, and 

75 the convolution with the filter are carried out on the two image signals S1 and S2. Thereafater, the 
subtraction processing is carried out. This operation is substantially equivalent to the operation, In which the 
subtraction is carried out such that the subtraction ratio may be different between the low frequency band 
and the high frequency band of the image signals. In this manner, the subtraction signal can be obtained 
such that the noise components obtained from the stimulable phosphor sheet, which was located at the 

20 position remote from the radiation source, may be reduced, and such that the information of the high 
frequency band obtained from the stimulable phosphor sheet, which was located at the position close to the 
radiation source, may be emphasized. Therefore, from the subtraction signal thus obtained, an Image 
having good image quality and containing little noise component can be reproduced. Also, the amount of 
calculation is not so large as in the wavelet transform or the Fourier transform. Therefore, the apparatus for 

25 carrying out the second embodiment of the energy subtraction processing method in accordance with the 
present invention can be kept simple. Further, the operation can be carried out quickly. 

A third embodiment of the superposition processing method for a radiation image in accordance with 
the present invention and a third embodiment of teh energy subtraction processing method in accordance 
with the present invention will be described herelnbelow. 

30 Image signals SI and S2 are obtained in the same manner as that described above with reference to 
Rgures 1 and 2. Thereafter, a superposition processing is carried out on the image signals SI and S2. 
Rgure 33 is a block diagram showing an apparatus for carrying out the third embodiment of the 
superposition processing method for a radiation image in accordance with the present invention. First, the 
image signals S1 and S2 are read from an image file 18A and an image file 18B in the storage medium 18 

35 and fed into an image processing means 219. The image processing with a filter described below is carried 
out on the two image signals S1 and S2. which have been fed into the image processing means 219. The 
image signals, which have been obtained from the image processing, are fed into an addition means 220 
and subjected to an addition. An addition signal, which has been obtained from the addition means 220. is 
fed into an image reproducing means 221, such as a CRT display device, and Is used for reproducing a 

40 visible image. 

How the image processing is carried out in the image processing means 219 will be described 
hereinbelow. 

As described above with reference to Figures 6A through 9, the optimum addition ratio, which gives the 
largest DQE serving as an index of an image, varies for different frequencies. Also, as described above, the 

45 image signal SI carries the information up to the high frequency band. However, in the image signal S2. 
noise components are predominant in the information of the high frequency band. Further, it is often desired 
that desired frequency components of the addition signal can be emphasized. Therefore, the convolution of 
the image signal SI and the convolution of the image signal S2 are carried out with filters, which are 
capable of emphasizing the response with respect to the high frequency band of the image signal Si, 

50 capable of reducing the response with respect to the high frequency band of the image signal S2, and 
capable of emphasizing the desired frequency components of the addition signal. Thereafter, the resulting 
image signals SI and S2 are added to each other. In this manner, an addition signal can be obtained which 
represents a radiation image having good image quality. 

How the Image signals 31 and 52 are processed will be described hereinbelow. 

55 As a filter for reducing the high frequency components of the image signal S2, the following filter F2 
may be employed. 
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The filter F2 has the frequency characteristics shown in Figure 34, 

Also, it may be considered to set a filter F1 so as to satisfy the condition of F1 + F2 = 1 . and the energy 
70 of the signal values of the addition signal may thereby be rendered identical with the energy of the original 
image signals Si and S2. Specifically, in such cases, the median (16/104) of the filter F2 is altered to a 
value of (1 - median), and teh signs of the other filter elements are reversed. In this manner, the filter F1 
having the values shown below may be obtained. 
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20 

The filter F1 has the frequency characteristics obtained by subtracting the characteristics shown in Figure 
34 from a value of 1 . 

In cases where the convolution of the image signal SI with the filter F1 and the convolution of the 
25 image signal S2 with the filter F2 may then be carried out, and an addition signal may be obtained. Also, 
the frequency characteristics of the addition signal, for example, at frequencies in the vicinity of 2 
cycles/mm, may then be emphasized. In such cases, it is necessary for the convolution of the addition 
signal to be carried out with a filter F3 having the frequency characteristics shown in Figure 35. In such 
cases, the amount of calculation becomes large, and a long time is required to obtain the addition signal. 
30 The values of the filter F3 are shown below. 
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Therefore, in the third embodiment of the superposition processing method for a radiation image in 
45 accordance with the present invention, the values of the filters F1 and F2 are altered such that the addition 
signal, which is obtained by adding the images signals resulting from the convolution of the image signal 31 
with the filter F1 and the convolution of the image signal S2 with the filter F2, may have the characteristics 
shown in Figure 35, i.e. such that F1 + F2 may have the frequency characteristics of the filter F3. 

Specifically, the convolution of the filter F1 with the filter F3 and the convolution of the filter F2 with the 
50 filter F3 are carried out with Formulas (28) and (29). 

Fr = Fi *F3 (28) 
Fa'^Fs-Fg (29) 

55 where FrF3 represents the convolution of the filter F1 with the fitter F3. In this manner, new filters FV and 
F2' are obtained. Thereafter, the convolution of the image signal S1 with the filter FV and the convolution of 
the image signal S2 with the filter F2' are carried out. 
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Specifically, the filter FV has the values shown below and the frequency characteristics shown in Figure 

36. 
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Also, the filter F2' has the values shown below and the frequency characteristics shown In Figure 37. 
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The filters FV and ¥Z are stored in the image processing nneans 219. The convolution of the image 
signal SI, which has been fed Into the image processing means 219. is earned out with the filter FV. Also, 
40 the convolution of the image signal S2, which has been fed into the image processing means 219, is carried 
out with the filter F2*. In this manner, processed image signals SV and S2' are obtained. The processed 
image signals SV and S2' are fed Into an addition means 220 and added to each other. 
The processing described above is represented by Formula (32). 

45 Sadd = Fi'*Si + F2''S2 (32) 

The addition signal Sadd obtained in this manner is fed into an image reproducing means 221 and used 
for reproducing a visible image. 

The image reproducing means 221 may be a display device, such as a cathode ray tube (CRT) display 
50 device, or a recording apparatus for recording an image on photosensitive film by scanning the photosen- 
sitive film with a light beam. 

As illustrated in Figure 38, in the visible image reproduced from the addition signal Sadd in the image 
reproducing means 221, the frequency characteristics in the vicinity of 2 cycles/mm have been em- 
phasized. 

55 In the manner described above, the convolution of the image signal SI is canied out with the filter FV 
described above. Also, the convolution of the image signal S2 is carried out with the filter F2' described 
above. The image signals obtained from the convolution are then added to each other. In this manner, the 
addition signal can be obtained such that the noise components obtained from the stimulable phosphor 
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sheet, which was located at the position remote from the radiation source, may be reduced, such that the 
information of the high frequency band obtained from the stimulable phosphor sheet, which was located at 
the position close to the radiation source, may be emphasized, and such that the addition signal may have 
the desired frequency characteristics. Therefore, from the addition signal thus obtained, an image having 

5 the desired frequency characteristics and good image quality and containing little noise component can be 
reproduced. Also, the amount of calculation is not so large as in the wavelet transform or the Fourier 
transform. Therefore, the apparatus for carrying out the third embodiment o1 the superposition processing 
method for a radiation image in accordance with the present invention can be kept simple. Further, the 
operation can be carried out quickly. 

70 In the aforesaid third embodiment of the superposition processing method for a radiation image in 
accordance with the present invention, the image processing Is canned out on both of the image signals S1 
and S2, Alternatively, the image processing described above may be carried out on either one of the image 
signals SI and S2. However, in cases where the image processing is carried out on both of the image 
signals SI and S2, an addition signal representing an image having better image quality can be obtained. 

75 Also, in the aforesaid third embodiment of the superposition processing method for a radiation image in 
accordance with the present invention, the image signals Si and 82 are obtained by scanning the 
stimulable phosphor sheets 4A and 4B, on which the radiation images have been stored, in the main 
scanning direction and sub-scanning direction. Therefore, the frequency characteristics of each of the image 
signals S1 and S2 vary for the main scanning direction and the sub-scanning direction. Specifically, with 

20 respect to the sub-scanning direction, in each of the image signals S1 and S2, the response becomes high 
for the frequency band due to the rotation speed of the scanning mirror 12 shown in Figure 2 and the 
speed, with which each of the stimulable phosphor sheets 4A and 4B is conveyed. 

Therefore, with respect to the sub-scanning direction, the convolution of each of the image signals 81 
and S2 is carried out with a filter capable of reducing the response for the aforesaid frequency band. Also, 

25 with respect to the main scanning direction, the convolution of each of the image signals SI and S2 is 
carried out with a filter capable of emphasizing the response for the aforesaid frequency band. The resulting 
image signals are then added to each other. From the addition signal thus obtained, an image can be 
reproduced which has the well-balanced frequency characteristics with respect to the main scanning 
direction and the sub-scanning direction, and which has good image quality and can serve as an effective 

30 tool in, particularly, the efficient and accurate diagnosis of an illness. 

In the aforesaid third embodiment of the superposition processing method for a radiation image in 
accordance with the present invention, in order for the image signals SI and S2 to be added to each other 
In a desired addition ratio, it is necessary for the value of the filter coefficient to be determined such that the 
frequency characteristics of the filter F2 may become the desired characteristics. The filter F2 may be 

35 determined with one of the two techniques for determining the filter F2, which have been described above 
with reference to the aforesaid modification of the second embodiment of the superposition processing 
method for a radiation image in accordance with the present invention. 

In the aforesaid third embodiment of the superposition processing method for a radiation Image in 
accordance with the present invention, as illustrated in Rgure 1, the radiation images are recorded on the 

40 two stimulable phosphor sheets 4A and 4B. The image signals obtained from the two stimulable phosphor 
sheets 4A and 48 are then added to each other. Alternatively, as illustrated in Figure 14, the radiation image 
of the object 1 may be recorded on the single stimulable phosphor sheet 4A. As illustrated in Figure 15. two 
Image signals to be added to each other may be obtained by detecting light emitted by the two surfaces of 
the stimulable phosphor sheet 4A. 

45 The two image signals SI and S2. which have been obtained from the analog-to-digital converters 17a 
and 17b in the same manner as that described above with reference to Figure 15, are added to each other 
in the same manner as that in the aforesaid third embodiment of the superposition processing method for a 
radiation image in accordance with the present invention. In this manner, as in the aforesaid third 
embodiment of the superposition processing method for a radiation image in accordance with the present 

50 invention, an image having the desired frequency characteristics and good image quality and containing 
little noise component can be reproduced from the addition signal, which has thus been obtained. 

As another alternative, as Illustrated in Figure 16, the two surfaces of the stimulable phosphor sheet 4A 
may be scanned respectively with laser beams 11a and lib, which have been produced by the laser beam 
sources 10a and 10b. The light emitted by the two surfaces of the stimulable phosphor sheet 4A may thus 

55 be detected, and two image signals may thereby be obtained. 

In the aforesaid third embodiment of the superposition processing method for a radiation image in 
accordance with the present invention, the two image signals are superposed one upon the other. The 
processing described above can also be carried out when energy subtraction processing is carried out on 
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two image signals. How the image processing is carried out on the image signals to be subjected to the 
energy subtraction processing will be described hereinbelow. 

In the image recording apparatus for carrying out the one-shot energy subtraction processing, which is 
shown in Rgure 20. In the same manner as that described above, the two radiation images, in which 
5 different images of a specific structure of the object 1 are embedded, are stored on the two stimulable 
phosphor sheets 4A and 4B. 

Thereafter, the radiation images are respectively read out from the two stimulable phosphor sheets 4A 
and 4B by using the radiation image read-out apparatus shown in Figure 2. Two digital image signals SI 
and S2 representing the radiation images are thereby obtained. The obtained image signals SI and S2 are 
70 stored in the storage medium 1 8. 

Subtraction processing is then carried out on the digital image signals S1 and S2, which have been 
obtained in the manner described above. Figure 39 Is a block diagram showing an apparatus for carrying 
out a third embodiment of the energy subtraction processing method in accordance with the present 
invention. First, the image signals Si and S2 are read from an image file 18A and an image file 18B in the 
15 storage medium 18 and fed into an image processing means 230. In the same manner as that in the 
aforesaid third embodiment of the superposition processing method for a radiation image In accordance 
with the present invention, the two image signals SI and S2 are subjected to the processing, such as the 
convolution with the mask, in the image processing means 230. Processed image signals SI' and S2* are 
thereby obtained. 

20 The image signals SV and S2' thus obtained are fed into a subtraction means 231. In the subtraction 
means 231, subtraction processing is carried out on the image signals Sr and S2\ 
Specifically, the subtractions are carried out with Formula (33) 

Ssub = ti • Fi * Si • - 12 • F2 • S2' (33) 

25 

where t1 and t2 represent the energy subtraction factors. 

In the manner described above, the subtraction signal Ssub Is obtained from the subtraction means 
231. The subtraction signal Ssub Is fed into an image reproducing means 232 and used for reproducing a 
visible image. 

30 In the aforesaid third embodiment of the energy subtraction processing method in accordance with the 
present invention, the one-shot energy subtraction processing is carried out. The aforesaid third embodi- 
ment of the energy subtraction processing method in accordance with the present invention is also 
applicable to the two-shot energy subtraction processing. 

In the manner described above, the convolution of the image signal SI is earned out with the filter F1. 

35 Also, the convolution of the image signal S2 is carried out with the filter F2. The image signals obtained 
from the convolution are then subtracted from each other. In this manner, the subtraction signal can be 
obtained such that the noise components obtained from the stimulable phosphor sheet, which was located 
at the position remote from the radiation source, may be reduced, such that the information of the high 
frequency band obtained from the stimulable phosphor sheet, which was located at the position dose to the 

40 radiation source, may be emphasized, and such that the subtraction signal may have the desired frequency 
characteristics. Therefore, from the subtraction signal thus obtained, an image having the desired frequency 
characteristics and good image quality and containing little noise component can be reproduced. Also, the 
amount of calculation is not so large as in the wavelet transform or the Fourier transfonm. Therefore, the 
apparatus for carrying out the third embodiment of the energy subtraction processing method In accordance 

45 with the present invention can be kept simple. Further, the operation can be carried out quickly. 

A fourth embodiment of the superposition processing method for a radiation image in accordance with 
the present invention and a fourth embodiment of the energy subtraction processing method in accordance 
with the present invention will be described hereinbelow. 

Image signals SI and S2 are obtained in the same manner as that described above with reference to 

50 Rgures 1 and 2. 

At this time, as illustrated in Figure 40, a radiation dose calculating means 319 calculates the dose of 
radiation delivered to the object 1 in accordance with the values of the output signals SA and SB and the 
read-out sensitivity and latitude in the photomultiplier 15. A signal SM representing the calculated dose of 
radiation is stored on the storage medium 18 together with the image signals SI and S2. 
55 Thereafter, a superposition processing is carried out on the Image signals Si and S2. Rgure 41 is a 
block diagram showing an apparatus for carrying out the fourth embodiment of the superposition processing 
method for a radiation image in accordance with the present invention. First, the image signals Si and S2 
and the signal SM are read from an Image file 18A and an image file 18B In the storage medium 18 and fed 
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into an image processing means 320. The image processing with a filter described below is carried out on 
the two image signals SI and S2, which have been fed Into the image processing means 320. The image 
signals, which have been obtained from the image processing, are fed into an addition means 321 and 
subjected to an addition. An addition signal Sadd, which has been obtained from the addition means 321. is 
5 fed Into an image reproducing means 322, such as a CRT display device, and Is used for reproducing a 
visible image. 

How the image processing is carried out in the image processing means 320 will be described 
hereinbelow. 

As described above with reference to Figures 6A through 9, the optimum addition ratio, which gives the 

10 largest DQE serving as an index of an image, varies for different frequencies. Also, as described above, the 
image signal S1 carries the information up to the high frequency band. However, in the image signal S2, 
noise components are predominant In the information of the high frequency band. Further, the frequency 
band, in which the noise components are predominant In the image signal, varies for different doses of 
radiation delivered to the object 1. Therefore, the dose of radiation delivered to the object 1 is taken into 

75 consideration, and a mask filter is obtained which is capable of emphasizing the response with respect to 
the high frequency band of the image signal SI. Also, a mask fitter is obtained which is capable of reducing 
the response with respect to the high frequency band of the image signal S2. The convolution of each of 
the image signals S1 and S2 is carried out with the corresponding mask filter. Thereafter, the resulting 
image signals SI and S2 are added to each other. In this manner, an addition signal can be obtained which 

20 represents a radiation Image having good Image quality. 

How the image signals Si and S2 are processed will be described hereinbelow. 
Rrst, the filters F1 and F2 to be used for the convolutions of the image signals Si and S2 are 
determined in accordance with the dose of radiation delivered to the object 1 . When the dose of radiation 
delivered to the object 1 is large, in the noise components, the proportion of the fixed noise, which is due to 

25 the structure of the stimulable phosphor sheet, such as the state in which the stimulable phosphor is 
applied to the stimulable phosphor sheet, becomes higher than the proportion of the quantum noise of the 
radiation. Specifically, between the image signals Si and S2, the correlation of the quantum noise of the 
radiation is small, but the correlation of the fixed noise is large. Therefore, when the dose of radiation 
delivered to the object 1 becomes large, the proportion of the noise components becomes high in the high 

30 frequency components of the image signal S2, which Is obtained from the stimulable phosphor sheet 4B 
located at the position remote from the radiation source during the image recording operation. Accordingly, 
as for the image signal S2, when the dose of radiation delivered to the object 1 is 0.1 mR. a filter F2 having 
the frequency characteristics shown In Figure 42 is set. Also, when the dose of radiation delivered to the 
object 1 is ImR, a filter F2 having the frequency characteristics shown in Figure 43 Is set. Further, when the 

35 dose of radiation delivered to the object 1 is lOmR, a filter F2 having the frequency characteristics shown in 
Rgure 44 is set. As for the Image signal SI. a filter F1 satisfying the condition of F1 + F2 = 1 is set such that 
the energy of the signal values of the addition signal may be identical with the energy of the original image 
signals SI and S2. The convolution of the image signal SI is carried out with the filter F1, and the 
convolution of the image signal S2 is carried out with the filter F2. 

40 The filter F2 having the frequency characteristics shown in Figure 42, the filter F2 having the frequency 
characteristics shown in Figure 43, and the filter F2 having the frequency characteristics shown in Figure 44 
can respectively be represented by Formulas (34), (35), and (36). 
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As for the filter F1 , such that it may satisfy the condition of F1 + F2 = 1 , the medians of the filters F2 
shown in Formulas (34), (35), and (36) are altered to a value of (1 - median), and the signs of the other filter 
elements are reversed. In this manner, the filters F1 are obtained which are represented by Fomnulas (37), 
(38), and (39). 
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lOgiolO- logioX 

F 1 (XmR) = . F I (1. OmR) + 

logiolO-logiol. 0 

logioX -logiol.O 

• F, (lOmR) 

logiolO-logiol. 0 



...(41) 

15 

In the same manner, the filter F2 is determined. 

After the filters F1 and F2 are detemnined in accordance with the dose of radiation delivered to the 
object 1 , the convolution of the image signal S1 is carried out with the filter F1 . Also, the convolution of the 
image signal S2 is carried out with the filter F2. In this manner, processed image signals SI' and S2' are 
20 obtained. The processed image signals SV and S2* are fed into an addition means 321 and added to each 
other. 

The processing described above is represented by Formula (42). 

Sadd = Fi * Si + F2 * S2 (42) 

25 

where FrS1 represents the convolution of SI with F1. 

The addition signal Sadd obtained in this manner is fed into an image reproducing means 322 and used 
for reproducing a visible image. 

The Image reproducing means 322 may be a display device, such as a cathode ray tube (CRT) display 

30 device, or a recording apparatus for recording an image on photosensitive film by scanning the photosen- 
sitive film with a light beam. 

In the manner described above, the convolution of the image signal Si is carried out with the filter F1 
described above. Also, the convolution of the image signal S2 is carried out with the filter F2 described 
above. The image signals obtained from the convolution are then added to each other. In this manner, the 

35 addition signal can be obtained such that the noise components obtained from the stimulable phosphor 
sheet, which was located at the position remote from the radiation source, may be reduced, such that the 
information of the high frequency band obtained from the stimulable phosphor sheet, which was located at 
the position close to the radiation source, may be emphasized, and such that the dose of radiation delivered 
to the object may be taken into consideration. Therefore, from the addition signal thus obtained, an image, 

40 which has good image quality and in which the noise components have been reduced in accordance with 
the dose of radiation delivered to the object, can be reproduced. Also, the amount of calculation is not so 
large as in the wavelet transform or the Fourier transform. Therefore, the apparatus for carrying out the 
fourth embodiment of the superposition processing method for a radiation image in accordance with the 
present invention can be kept simple. Further, the operation can be carried out quickly. 

45 In the aforesaid fourth embodiment of the superposition processing method for a radiation image in 
accordance with the present invention, the image processing is carried out on both of the image signals S1 
and S2. Alternatively, the image processing described above may be carried out on either one of the image 
signals SI and S2. However, in cases where the image processing is carried out on both of the image 
signals SI and S2, an addition signal representing an image having better image quality can be obtained. 

50 For reasons of the apparatus, for example, in cases where the capacity of the memory in the apparatus 
is small, it may be desired to carry out the convolution of the image signal S1 and the convolution of the 
image signal S2 with a single filter. In such cases, the addition may be carried out with the formulas shown 
below. 

55 Sadd = (Si - F2 • Si) + (Fa ' S2) (43) 

(Because F1 + F2 = 1 , , 
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Sadd = F, • Si + F2 * S2 
= (1 - F2) ■ Si + Fz * S2 
= (S^ -F2*Si) + (F2*S2)) 
or Sadd = F, ' Si + S2 - Fi ' S2 (44) 

In this manner, the same resutts as those with Formula (42) can be obtained. 

In the aforesaid fourth embodiment of the superposition processing method for a radiation image in 
accordance with the present invention, in order for the Image signals SI and S2 to be added to each other 
in a desired addition ratio, it is necessary for the value of the filter coefficient to be determined such that the 
frequency characteristics of the filter F2 may become the desired characteristics. The filter F2 may be 
determined with one of the two techniques for determining the filter F2, which have been described above 
with reference to the aforesaid modification of the second embodiment of the superposrtion processing 
method for a radiation image in accordance with the present Invention. 

Also, in the aforesaid fourth embodiment of the superposition processing method for a radiation image 
in accordance with the present invention, the entire dose of radiation delivered to the object 1 is calculated, 
and the filters FI and F2 are set in accordance with the calculated dose of radiation. Alternatively, the dose 
of radiation delivered to the object 1 may be calculated for each of different portions of the object 1 , the 
patterns of which are embedded in the radiation image represented by each of the image signals S1 and 
S2. in accordance with, for example, the values of the respective image signal components of each of the 
image signals S1 and S2 or the unsharp mask signal of each of the image signals S1 and S2. Each of the 
filters FI and F2 may be set for each of the portions of the object in accordance with the calculated dose of 
radiation. Thereafter, the convolution of each of the image signals S1 and S2 with the filter, which has been 
set for each of the portions of the object. 

In cases where each of the filters F1 and F2 is set for each of portions of the object, the patterns of 
which are embedded in the radiation image, and the convolutions and the addition are then carried out, an 
addition signal can be obtained in an addition ratio optimum for each of the portions of the object. 
Therefore, it becomes possible to obtain the addition signal representing a radiation image which has good 
image quality with respect to each of the portions of the object and can serve as an effective tool in, 
particularly, the efficient and accurate diagnosis of an illness at each of the portions of the object. 

Further, in the aforesaid fourth embodiment of the superposition processing method for a radiation 
image in accordance with the present invention, the filters satisfying the condition of F1+F2 = 1 are 
employed. Alternatively, filters, which do not satisfy the condition of FI + F2 = 1 , may be employed. 
However, in cases where such filters are used, it is necessary to carry out a predetermined weighting 
process such that the energy of the signal values of the addition signal Sadd may become identical with the 
energy of the original image signals SI and S2. Specifically, it is necessary for the addition signal Sadd to 
be calculated with Formula (45). 

Sadd = t ♦ Fi * Si + (1 - 1) • F2 * S2 (45) 

In the aforesaid fourth embodiment of the superposition processing method for a radiation image in 
accordance with the present invention, as illustrated in Rgure 1. the radiation images are recorded on the 
two stimulable phosphor sheets 4A and 4B. The Image signals obtained from the two stimulable phosphor 
sheets 4A and 4B are then added to each other. Alternatively, as illustrated in Figure 14, the radiation image 
of the object 1 may be recorded on the single stimulable phosphor sheet 4A. As illustrated in Figure 15, two 
image signals to be added to each other may be obtained by detecting light emitted by the two surfaces of 
the stimulable phosphor sheet 4A. 

The two image signals SI and S2, which have been obtained from the analog-to-digital converters 17a 
and 17b in the same manner as that described above with reference to Figure 15, are added to each other 
in the same manner as that in the aforesaid fourth embodiment of the superposition processing method for 
a radiation image in accordance with the present invention. In this manner, as in the aforesaid fourth 
embodiment of the superposition processing method for a radiation image in accordance with the present 
invention, an image having good image quality and containing little noise component can be reproduced 
from the addition signal, which has thus been obtained. 

As another alternative, as illustrated in Figure 16. the two surfaces of the stimulable phosphor sheet 4A 
may be scanned respectively with laser beams 11a and lib, which have been produced by the laser beam 
sources 10a and 10b. The light emitted by the two surfaces of the stimulable phosphor sheet 4A may thus 
be detected, and two image signals may thereby be obtained. 
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In the aforesaid fourth embodiment of the superposition processing method for a radiation image in 
accordance with the present invention, the two image signals are superposed one upon the other. The 
processing described above can also be carried out when energy subtraction processing is carried out on 
two image signals. How the image processing is carried out on the image signals to be subjected to the 

5 energy subtraction processing will be described hereinbelow. 

In the image recording apparatus for carrying out the one-shot energy subtraction processing, which is 
shown in Rgure 20. in the same manner as that described above, the two radiation images, in which 
different images of a specific structure of the object 1 are embedded, are stored on the two stimulable 
phosphor sheets 4A and 4B. 

70 Thereafter, the radiation images are respectively read out from the two stimulable phosphor sheets 4A 
and 4B by using the radiation image read-out apparatus shown in Figure 2. Two digital image signals S1 
and S2 representing the radiation images are thereby obtained. The obtained image signals S1 and S2 are 
stored in the storage medium 18. 

Subtraction processing is then carried out on the digital image signals S1 and S2, which have been 

75 obtained in the manner described above. Figure 45 is a block diagram showing an apparatus for carrying 
out a fourth embodiment of the energy subtraction processing method in accordance with the present 
invention. First, the image signals SI and S2 are read from an image file 18A and an image file 18B in the 
storage medium 18 and fed into an image processing means 330. In the same manner as that in the 
aforesaid fourth embodiment of the superposition processing method for a radiation image in accordance 

20 with the present invention, the two image signals SI and S2 are subjected to the processing, such as the 
convolution with the mask, in the image processing means 330. Processed image signals SV and S2* are 
thereby obtained. 

The image signals SI * and S2' thus obtained are fed into a subtraction means 331 . In the subtraction 
means 331. subtraction processing is carried out on the image signals SV and S2\ 
25 Specifically, the subtractions are carried out with Formula (46) 

Ssub = ti • Fi *Si'-t2 • Fz •S2' (46) 

where t1 and t2 represent the energy subtraction factors. 
30 In the manner described above, the subtraction signal Ssub is obtained from the subtraction means 
331 . The subtraction signal Ssub is fed into an image reproducing means 332 and used for reproducing a 
visible image. 

In the aforesaid fourth embodiment of the energy subtraction processing method in accordance with the 
present invention, the one-shot energy subtraction processing is carried out. The aforesaid fourth embodi- 

35 ment of the energy subtraction processing method in accordance with the present invention is also 
applicable to the two-shot energy subtraction processing. 

In the manner described above, the convolutions of the image signals SI and S2 are carried out with 
the filters described above. The image signals obtained from the convolution are then subtracted from each 
other. In this manner, the difference signal can be obtained such that the noise components obtained from 

40 the stimulable phosphor sheet, which was located at the position remote from the radiation source, may be 
reduced, such that the information of the high frequency band obtained from the stimulable phosphor sheet, 
which was located at the position close to the radiation source, may be emphasized, and such that the dose 
of radiation delivered to the object may be taken into consideration. Therefore, from the difference signal 
thus obtained, an image, which has good image quality and in which the noise components have been 

45 reduced in accordance with the dose of radiation delivered to the object, can be reproduced. Also, the 
amount of calculation is not so large as in the wavelet transform or the Fourier transform. Therefore, the 
apparatus for carrying out the fourth embodiment of the energy subtraction processing method in 
accordance with the present invention can be kept simple. Further, the operation can be carried out quickly. 

50 Claims 

1. A superposition processing method for a radiation Image, comprising the steps of: 

i) obtaining a plurality of image signals, which represent a radiation image of a single object or 
radiation images of the single object, and which have different frequency characteristics, each of the 

55 image signals being made up of a series of image signal components, 

ii) weighting the image signal components of the plurality of the image signals, which image signal 
components represent conresponding picture elements, with predetermined weight factors, and 
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iii) adding the weighted image signal cxjmponents of the plurality of the innage signals to one 
another, which image signal components represent corresponding picture elements, an addition 
signal being thereby obtained, 

wherein the value of the weight factor with respect to the frequency components, which have a 
6 low signai-to-noise ratio, is rendered smaller than the value of the weight factor with respect to the 

frequency components, which have a high signal-to-noise ratio, in accordance with the frequency 
characteristics of each of the image signals. 

2. A superposition processing method for a radiation image as defined in Claim 1 wherein the image 
70 superposition processing is carried out by: 

a) can7ing out a Fourier transfonm on each of the image signals, each of the image signals being 
thereby decomposed into a plurality of Fourier transfonm factor signals, each being of one of 
different frequency bands, 

b) weighting the Fourier transform factor signals, which are of a single frequency band, with the 
75 weight factors, such that the values of the weight factors may be varied for the Fourier transfonm 

factor signals of the different frequency bands, 

c) adding the weighted Fourier transform factor signals, which are of a single frequency band, to one 
another, an addition Fourier transform factor signal being thereby obtained for each of the different 
frequency bands, and 

20 d) carrying out an inverse Fourier transform on the addition Fourier transform factor signal. 

3. A superposition processing method for a radiation image as defined in Claim 1 wherein the image 
superposition processing is carried out by: 

a) subjecting each of tiie image signals to a transform to multi-resolution space, each of the image 
25 signals being thereby decomposed into a plurality of transform factor signals, each being of one of 

different frequency bands, 

b) weighting the transfomn factor signals, which are of a single frequency band, with the weight 
factors, such that the values of the weight factors may be varied for the transform factor signals of 
the different frequency bands, 

30 c) adding the weighted transform factor signals, which are of a single frequency band, to one 

another, an addition transform factor signal being thereby obtained for each of the different 
frequency bands, and 

d) carrying out an inverse ti-ansform on the addition transform factor signal. 

35 4. A superposition processing metfiod for a radiation image as defined in Claim 3 wherein the transform to 
multi-resolution space is carried out with a wavelet transform. 

5. A superposition processing method for a radiation image as defined in Claim 1, 3, or 4 wherein the 
value of the weight factor with respect to a portion of the radiation image, to which a large dose of 
40 radiation reached during an operation for recording the radiation image, is rendered larger than the 
value, of the weight factor with respect to a portion of the radiation image, to which a small dose of 
radiation reached during the operation for recording the radiation image, in accordance with the 
portions of the object, the patterns of which are embedded in tiie radiation image. 

45 6. A superposition processing method for a radiation image as defined in Claim 1 wherein the plurality of 
the image signals are analog image signals, and the image superposition processing is carried out by: 

a) filtering all of the analog image signals with filters, which have weights changing the frequency 
characteristics of the analog image signals, and 

b) adding the analog image signals, v/hich have been obtained from the filtering process, to one 
50 another. 

7. A superposition processing method for a radiation image as defined in Claim 1 wherein the plurality of 
the image signals are analog image signals, and the image superposition processing is canned out by: 

a) tiltering at least a single desired analog image signal, which is among the analog image signals. 
55 with a filter, which has a weight changing the frequency characteristics of the desired analog image 

signal, and 

b) adding the analog image signal, which has been obtained from the tiltering process, and the other 
analog image signals to one another. 
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8. A superposition processing method for a radiation image as defined in Claim 1 , 2. 3, 4, 6, or 7 wherein 
the plurality of the image signals are obtained by: 

exposing the two surfaces or either one of the two surfaces of a single stimulable phosphor sheet, 
on which the radiation image has been stored, to stimulating rays, which cause the stimulable phosphor 
sheet to emit light in proportion to the amount of energy stored thereon during its exposure to radiation, 
and 

photoelectrically detecting the emitted light independently on the opposite surface sides of the 
stimulable phosphor sheet. 

9. A superposition processing method for a radiation image as defined in Claim 5 wherein the plurality of 
the image signals are obtained by: 

exposing the two surfaces or either one of the two surfaces of a single stimulable phosphor sheet, 
on which the radiation image has been stored, to stimulating rays, which cause the stimulable phosphor 
sheet to emit light in proportion to the amount of energy stored thereon during its exposure to radiation, 
and 

photoelectrically detecting the emitted light independently on the opposite surface sides of the 
stimulable phosphor sheet. 

10. A superposition processing method for a radiation image as defined in Claim 1, 2, 3, 4, 6, or 7 wherein 
the plurality of the image signals are obtained by: 

exposing each of at least two stimulable phosphor sheets, on each of which a radiation image has 
been stored, to stimulating rays, which cause the stimulable phosphor sheet to emit light in proportion 
to the amount of energy stored thereon during its exposure to radiation, and 

photoelectrically detecting the emitted light. 

11. A superposition processing method for a radiation image as defined in Claim 5 wherein tiie plurality of 
the image signals are obtained by: 

exposing each of at least two stimulable phosphor sheets, on each of which a radiation image has 
been stored, to stimulating rays, which cause the stimulable phosphor sheet to emit light in proportion 
to the amount of energy stored thereon during its exposure to radiation, and 

photoelectrically detecting the emitted light. 

12. An energy subtraction processing method comprising the steps of: 

i) forming a plurality of radiation images of a single object respectively with a plurality of kinds of 
radiation having different energy distributions, different images of at least part of the object being 
embedded in the plurality of the radiation images, 

ii) obtaining a plurality of image signals, each of which is made up of a series of image signal 
components, from the plurality of the radiation images, 

iii) weighting the image signal components of the plurality of the image signals, which image signal 
components represent corresponding picture elements in the plurality of the radiation images, with 
predetermined weight factors, and 

iv) subtracting the weighted image signal components of the plurality of the image signals from one 
another, which image signal components represent corresponding picture elements in the plurality of 
the radiation images, a difference signal representing an image of a specific structure of the object 
being thereby obtained. 

wherein the value of the weight factor with respect to the frequency components, which have a 
low signal-to-noise ratio, is rendered smaller than the value of the weight factor with respect to the 
frequency components, which have a high signal-to-noise ratio, in accordance with the frequency 
characteristics of each of the image signals. 

13. An energy subtraction processing method as defined in Claim 12 wherein the difference signal is 
obtained by: 

a) carrying out a Fourier transform on each of the image signals, each of the image signals being 
thereby decomposed into a plurality of Fourier transform factor signals, each being of one of 
different frequency bands, 

b) weighting the Fourier transform factor signals, which are of a single frequency band, with the 
weight factors, such that tine values of the weight factors may be varied for the Fourier transform 
factor signals of the different frequency bands, 
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c) subtracting the weighted Fourier transform factor signals, which are of a single frequency band, 
from one another, a subtraction Fourier transform factor signal being thereby obtained for each of 
the different frequency bands, and 

d) carrying out an inverse Fourier transform on the subtraction Fourier transform factor signal, 

5 

14. An energy subtraction processing method as defined in Claim 12 wherein the subtraction signal is 
obtained by: 

a) subjecting each of the image signals to a transform to multi-resolution space, each of the image 
signals being thereby decomposed into a plurality of transform factor signals, each being of one of 

70 different frequency bands, 

b) weighting the transfonm factor signals, which are of a single frequency band, with the weight 
factors, such that the values of the weight factors may be varied for the transfonm factor signals of 
the different frequency bands, 

c) subtracting the weighted transform factor signals, which are of a single frequency band, from one 
;5 another, a subtraction transform factor signal being thereby obtained for each of the different 

frequency bands, and 

d) carrying out an inverse transform on the subtraction transform factor signal. 

15. An energy subtraction processing method as defined in Claim 14 wherein the transform to multi- 
20 resolution space is carried out with a wavelet transform. 

16. An energy subtraction processing method as defined in Claim 12, 14, or 15 wherein the value of the 
weight factor with respect to a portion of each radiation image, to which a large dose of radiation 
reached during an operation for recording the radiation image, is rendered larger than the value of the 

25 weight factor with respect to a portion of the radiation image, to which a small dose of radiation reached 
during the operation for recording the radiation image, in accordance with the portions of the object, the 
patterns of which are embedded in the radiation image. 

17. An energy subtraction processing method as defined In Claim 12 wherein the plurality of the image 
signals are analog image signals, and the difference signal is obtained by: 

a) filtering all of the analog Image signals with filters, which have weights changing the frequency 
characteristics of the analog image signals, and 

b) subtracting the analog image signals, which have been obtained from the filtering process, from 
one another. 

18. An energy subtraction processing method as defined in Claim 12 wherein the plurality of the image 
signals are analog image signals, and the difference signal is obtained by: 

a) filtering at least a single desired analog image signal, which is among the analog image signals, 
with a filter, which has a weight changing the frequency characteristics of the desired analog image 
signal, and 

b) subtracting the analog image signal, which has been obtained from the filtering process, and the 
other analog image signals from one another. 

19. A superposition processing method for a radiation image, comprising the steps of: 

i) obtaining a plurality of image signals, which represent a radiation image of a single object or 
radiation images of the single object, and which have different frequency characteristics, each of the 
image signals being made up of a series of image signal components, and 

ii) adding the image signal components of the plurality of the image signals to one another, which 
image signal components represent corresponding picture elements, an addition signal being 
thereby obtained, 

wherein, such that the signal-to-noise ratio of the addition signal may be kept high, the addition 
signal is obtained by: 

a) subjecting at least a single desired image signal, which is among the plurality of the image 
signals, to image processing, which changes the frequency characteristics of the desired image 
signal, and 

b) adding the image signal, which has been obtained from the image processing, and the other 
Image signals to one another. 
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20. A superposition processing method for a radiation image as defined in Claim 19 wherein the image 
processing is the processing for carrying out the convolution of the desired image signal with a mask 
filter, which has predetermined frequency characteristics. 

5 21. A superposition processing method for a radiation image as defined in Claim 20 wherein the sum of the 
frequency characteristics of the mask filters employed for the plurality of the Image signals is equal to 
1 at an arbitrary frequency. 

22. A superposition processing method for a radiation Image as defined in Claim 19 wherein each of the 
70 plurality of the image signals is obtained by scanning a stimulable phosphor sheet, on which the 

radiation image has been stored, with stimulating rays, which cause the stimulable phosphor sheet to 
emit light in proportion to the amount of energy stored thereon during its exposure to radiation, in the 
main scanning direction and the sub-scanning direction, the emitted light being detected photoelec- 
trically, and 

75 the image processing Is the processing for emphasizing the high frequency components of the 

desired image signal with respect to the main scanning direction. 

23. A superposition processing method for a radiation image as defined in Claim 22 wherein said 
processing is the processing for carrying out the convolution of the desired image signal with respect to 

20 the main scanning direction with a mask filter, which has frequency characteristics capable of 
emphasizing the high frequency components of the desired image signal with respect to the main 
scanning direction. 

24. A superposition processing method for a radiation image as defined in Claim 19, 20, 21, 22, or 23 
25 wherein all of the plurality of the image signals are subjected to the image processing, which changes 

the frequency characteristics. 

25. A superposition processing method for a radiation image, comprising the steps of: 

i) obtaining two image signals, which represent a radiation image of a single object or radiation 
30 images of the single object, and which have different frequency characteristics, each of the image 

signals being made up of a series of image signal components, and 

ii) adding the image signal components of the two image signals to one another, which image signal 
components represent corresponding picture elements, an addition signal being thereby obtained, 

wherein, such that the signal-to-noise ratio of the addition signal may be kept high, the addition 
35 signal is obtained by: 

a) obtaining a difference signal between the two image signals, 

b) canrying out the convolution of the difference signal with either one of two mask filters, which have 
predetermined frequency characteristics, the sum of the frequency characteristics of the two mask 
filters being equal to 1 at an arbitrary frequency, and 

40 c) adding the difference signal, which has been obtained from the convolution, and either one of the 

two image signals to each other. 

26. An energy subtraction processing method comprising the steps of: 

j) forming a plurality of radiation Images of a single object respectively with a plurality of kinds of 
45 radiation having different energy distributions, different images of at least part of the object being 

embedded in the plurality of the radiation images, 

Ii) obtaining a plurality of image signals, each of which is made up of a series of image signal 
components, from the plurality of the radiation images, 

Hi) weighting the image signal components of the plurality of the image signals, which image signal 
50 components represent corresponding picture elements in the plurality of the radiation images, with 

predetermined weight factors, and 

Iv) subtracting the weighted image signal components of the plurality of the image signals from one 
another, which image signal components represent conresponding picture elements in the plurality of 
the radiation images, a difference signal representing an image of a specific structure of the object 
55 being thereby obtained, 

wherein, such that the signal-to-noise ratio of the difference signal may be kept high, the 
difference signal Is obtained by: 
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a) subjecting at least a single desired image signal, which is among the plurality of the Image 
signals, to image processing, which changes the frequency characteristics of the desired image 
signal, and 

b) subtracting the image signal, which has been obtained from the image processing, and the other 
5 image signals from one another. 

27. An energy subtraction processing method as defined in Claim 26 wherein the image processing is the 
processing for carrying out the convolution of the desired image signal with a predetermined mask 
filter. 

28. An energy subtraction processing method as defined in Claim 26 or 27 wherein all of the plurality of the 
image signals are subjected to the image processing, which changes the frequency characteristics. 

29. A superposition processing method for a radiation image, comprising the steps of: 

/5 i) obtaining a plurality of image signals, which represent a radiation image of a single object or 

radiation images of the single object, and which have different frequency characteristics, each of the 
image signals being made up of a series of image signal components, and 

it) adding the image signal components of the plurality of the image signals to one another, which 
image signal components represent corresponding picture elements, an addition signal being 
20 thereby obtained, 

wherein the addition signal Is obtained by: 

a) carrying out the convolution of at least a single desired image signal, which is among the plurality 
of the image signals, with a mask filter, which has frequency characteristics capable of keeping the 
signal-to-noise ratio of the addition signal high and altering the response characteristics of the 

25 addition signal with respect to a desired frequency band when the addition signal is obtained by 

adding an image signal, that is obtained from the convolution of the desired Image signal, and the 
other image signals to one another, a processed image signal being thereby obtained from the 
convolution, and 

b) adding the processed image signal and the other image signals to one another, 

30 

30. A superposition processing method for a radiation image as defined in Claim 29 wherein each of the 
plurality of the image signals is obtained by scanning a sheet-like recording medium, on which the 
radiation image has been recorded, with a light beam in the main scanning direction and the sub- 
scanning direction, the radiation image being thereby read out photoelectrically. and 

35 the processed image signal is obtained by carrying out the convolution of the desired image signal 

with the mask filter, which has different frequency characteristics with respect to the main scanning 
direction and the sub-scanning direction, such that the frequency response characteristics of the 
addition signal with respect to the main scanning direction and the frequency response characteristics 
of the addition signal with respect to the sutj-scanning direction may become approximately identical 

40 with each other. 

31. A superposition processing method for a radiation image as defined in Claim 29 or 30 wherein the 
processing with the convolution is carried out on all of the plurality of the image signals. 

45 32. An energy subtraction processing method comprising the steps of: 

i) forming a plurality of radiation images of a single object respectively with a plurality of kinds of 
radiation having different energy distributions, different images of at least part of the object being 
embedded in the plurality of the radiation images, 

ii) obtaining a plurality of image signals, each of which is made up of a series of image signal 
50 components, from the plurality of the radiation images, 

ili) weighting the image signal components of the plurality of the image signals, which image signal 
components represent corresponding picture elements in the plurality of the radiation images, with 
predetermined weight factors, and 

Iv) subtracting the weighted image signal components of the plurality of the image signals from one 
55 another, which image signal components represent conresponding picture elements in the plurality of 

the radiation images, a difference signal representing an image of a specific structure of the object 
being thereby obtained. 

wherein the difference signal is obtained by: 
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a) carrying out the convolution of at least a single desired innage signal, which is annong the plurality 
of the Image signals, with a nnask filter, which has frequency characteristics capable of keeping the 
signal-to-noise ratio of the difference signal high and altering the response characteristics of the 
difference signal with respect to a desired frequency band when the difference signal is obtained by 

5 subtracting an image signal, that is obtained from the convolution of the desired image signal, and 

the other Image signals from one another, a processed image signal being thereby obtained from 
the convolution, and 

b) subtracting the processed image signal and the other Image signals from one another. 

70 33. An energy subtraction processing method as defined in Claim 32 wherein each of the plurality of the 
image signals is obtained by scanning a sheet-like recording medium, on which the radiation image has 
been recorded, with a light beam in the main scanning direction and the sub-scanning direction, the 
radiation image being thereby read out photoelectrically, and 

the processed image signal is obtained by carrying out the convolution of the desired image signal 

15 with the mask filter, which has different frequency characteristics with respect to the main scanning 
direction and the sub-scanning direction, such that the frequency response characteristics of the 
difference signal with respect to the main scanning direction and the frequency response characteristics 
of the difference signal with respect to the sub-scanning direction may become approximately identical 
with each other. 

20 

34. An energy subtraction processing method as defined in Claim 32 or 33 wherein the processing with the 
convolution is carried out on all of the plurality of the image signals, 

35. A superposition processing method for a radiation image, comprising the steps of: 

25 i) obtaining a plurality of image signals, which represent a radiation image of a single object or 

radiation images of the single object having been formed by delivering radiation to the single object, 
and which have different frequency characteristics, each of the image signals being made up of a 
series of Image signal components, and 

ii) adding the Image signal components of the plurality of the image signals to one another, which 
30 image signal components represent corresponding picture elements, an addition signal being 

thereby obtained, 

wherein the addition signal is obtained by: 

a) calculating the dose of radiation delivered to the object, 

b) setting a mask filter for at least a single desired image signal, which Is among the plurality of the 
35 image signals, in accordance with the calculated dose of radiation, the mask filter having frequency 

characteristics capable of keeping the signal-to-noise ratio of the addition signal high when the 
addition signal is obtained by adding an image signal, that is obtained from the convolution of the 
desired image signal, and the other image signals to one another, 

c) carrying out the convolution of the desired image signal with the mask filter, and 

40 d) adding the image signal, which has been obtained from the convolution, and the other image 

signals to one another. 

36. A superposition processing method for a radiation image as defined in Claim 35 wherein the sum of the 
frequency characteristics of the mask filters, which are employed for the plurality of the image signals, 

45 is equal to 1 at an arbitrary frequency. 

37. A superposition processing method for a radiation Image as defined in Claim 35 or 36 wherein the dose 
of radiation is calculated for each of portions of the object, the patterns of which are embedded in the 
radiation image. 

50 the mask filter is set for each of the portions of the object in accordance with the calculated dose of 

radiation, and 

the convolution of the desired image signal is carried out with the mask filter, which has been set 
for each of the portions of the object. 

55 38, A superposition processing method for a radiation image as defined in Claim 35 or 36 wherein the 
processing with the convolution is canried out on all of the plurality of the image signals. 
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39. A superposition processing method for a radiation image as defined in Claim 37 wherein the processing 
with the convolution is canried out on all of the plurality of the image signals. 

40. An energy subtraction processing method comprising the steps of: 

i) forming a plurality of radiation images of a single object respectively with a plurality of kinds of 
radiation having different energy distributions, different images of at least part of the object being 
embedded in the plurality of the radiation Images, 

ii) obtaining a plurality of image signals, each of which is made up of a series of image signal 
components, from the plurality of the radiation images. 

iil) weighting the image signal components of the plurality of the image signals, which image signal 
components represent corresponding picture elements in the plurality of the radiation images, with 
predetermined weight factors, and 

iv) subtracting the weighted image signal components of the plurality of the image signals from one 
another, which image signal components represent coaesponding picture elements in the plurality of 
the radiation images, a difference signal representing an image of a specific structure of the object 
being thereby obtained, 

wherein the difference signal is obtained by: 

a) calculating the dose of radiation delivered to the object. 

b) setting a mask filter for at least a single desired image signal, which is among the plurality of the 
image signals, in accordance with the calculated dose of radiation, the mask filter having frequency 
characteristics capable of keeping the signal-to-noise ratio of the difference signal high when the 
difference signal Is obtained by subtracting an image signal, that is obtained from the convolution of 
the desired image signal, and the other image signals from one another, 

c) carrying out the convolution of the desired image signal with the mask filter, and 

d) subtracting the image signal, which has been obtained from the convolution, and the other image 
signals from one another. 

41. An energy subtraction processing method as defined in Claim 40 wherein the dose of radiation is 
calculated for each of portions of the object, the patterns of which are embedded In the radiation image, 

the mask filter is set for each of the portions of the object in accordance with the calculated dose of 
radiation, and 

the convolution of the desired image signal is carried out with the mask filter, which has been set 
for each of the portions of the object. 

42. An energy subtraction processing method as defined in Claim 40 or 41 wherein the processing with the 
convolution is carried out on all of the plurality of the image signals. 

43. An energy subtraction processing method as defined in Claim 42 wherein the sum of the frequency 
characteristics of the mask filters, which are employed for the plurality of the image signals, is equal to 
1 at an arbitrary frequency. 
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